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Abstract. The existence of periodic waves propagating downstream on the 
surface of a two-dimensional infinitely deep water under gravity is established 
for a general class of vorticities. When reformulated as an elliptic boundary 
value problem in a fixed semi-infinite strip with a parameter, the operator 
describing the problem is nonlinear and non-Fredholm. A global connected 
set of nontrivial solutions is obtained via singular theory of bifurcation. Each 
solution on the continuum has a symmetric and monotone wave profile. The 
proof uses a generalized degree theory, global bifurcation theory and Wyburn's 
lemma in topology, combined with the Schauder theory for elliptic problems 
and the maximum principle. 

1. Introduction 

The problem of surface water waves, in its simplest form, concerns the two- 
dimensional dynamics of an incompressible inviscid fluid of infinite depth and the 
wave motion on its surface layer, acted upon by gravity. The effect of surface ten- 
sion is neglected. Suppose for definiteness that in the (x, y)-Cartesian coordinates 
gravity acts in the negative y-direction and that the fluid at time t occupies the 
region bounded above by the moving surface, given as the graph y = rj(t, x). In the 
fluid region the velocity field (u(t, x, y), v(t, x, y)) and the pressure P(t, x, y) satisfy 
the Euler equations 

U X +Vy = 0, U t + UU X + VUy = -P X , 

(1-1) 

v t + uv x + vv y = —P v - g, 

where g > denotes the acceleration due to gravity. Throughout, subscripts de- 
note partial derivatives. The flow is allowed to be rotational, characterized by the 
vorticity u> = v x — u y . The kinematic and dynamic boundary conditions 

(1.2) v = T]t + ur] x and P = P atm 

at the surface layer express, respectively, that the surface moves with the velocity 
of the fluid particles at the surface and that the pressure at the surface is constant 
atmospheric, denoted by P a tm- The boundary condition at the infinite bottom 

(1.3) (it, v) — * (0, 0) as y — * — oo 

states that the flow at great depths is practically at rest. 

It is a matter of common experience that waves which may be typically observed 
on the surface of the ocean or the river are approximately periodic and propagating 
of permanent form at a constant speed. Waves of this kind are referred to as Stokes 
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waves. Intuitively, they are symmetric waves whose profile rises and falls exactly 
once per wavelength. 

Waves of Stokes' kind are among the few exact solutions of the water-wave prob- 
lenfl l[LT ]) -([0 ]) . To clarify their existence is therefore a fundamental mathematical 
issue, and it is the subject of investigation here. Furthermore, they are a genuine 
nonlinear phenomenon. In his formal yet far-reaching consideration [Sto80] , Stokes 
observed that characteristics of actual water waves deviate significantly from what 
the linear theory predicts. In particular, he conjectured that a periodic wave of 
maximum height exists and it is distinguished by a sharp cusp at the wave crest 
with the contained angle of 120 degreqj- 

When uj = 0, namely, in the irrotational setting, the existence theory of Stokes 
waves dates back to the construction of small-amplitude waves by Levi-Civita 
[LC25] and independently by Nekrasov [Nek51], and it includes the global the- 
ory by Krasovskii [Kra61] and by Keady and Norbury [KN78]. Stokes' conjec- 
ture is proved in the works of Amick, Fracnkcl and Toland [AFT82, Tol78] and 
McLcod [McL79] . These results arc based on the reformulation of the problem as 
Nckrasov's integral equation and they are reviewed in [Tol96]. Further advances 
are made in [BDTOOa, BDTOOb, BT03] based on the formulation of the problem as 
Babenko's pseudo-differential equation [Bab87]. 

While the zero-vorticity setting may serve as an approximation under certain 
circumstances and it dominates in the existing literature, ocean currents typically 
carry vorticitics in their wind-drifted boundary layers [Mei89]. Moreover, the gov- 
erning equations of water waves allow for rotational motions. Gerstner [Ger09] in 
1809 (earlier than Stokes!) found an explicit formula for periodic traveling waves 
on deep water with a particular nonzero vorticity; in the irrotational setting, no 
traveling- wave solution of the water-wave problem (|1.1|) - (|1.3|) is known in the closed 
form. Perhaps, more striking are periodic traveling waves over a current with a crit- 
ical layer and with a closed streamline [Wah09, CV09]; such waves cannot exist in 
the irrotational setting. Vorticity, in addition, has subtle influence on the hydro- 
dynamic stability [HL08] of traveling water waves. 

The existence theory of irrotational Stokes waves is so quite complete that it is 
featured in a textbook [BT03], but their siblings, with vorticity, enjoy a different 
reputation. Dubreil-Jacotin [D J34] used a partial hodograph transform to reformu- 
late the original free boundary problem in a fixed domain and she addressed the 
existence of small-amplitude waves by power series methods. Zcidler [Zci73] later 
suggested to use a quasi-conformal mapping and treated small-amplitude waves 
in the finite-depth case as well as under gravity and surface tension. The global 
theory began only recently when Constantin and Strauss [CS04] recognized that 
Dubrcil-Jacotin's formulation could be regarded as an abstract operator equation 
in a Banach space. They employed the topological degree theory, as adapted by 
Healey and Simpson [HS98] for a general class of nonlinear elliptic operators, and 
the global bifurcation theory of Rabinowitz [Rab71], and in the finite-depth case 
they obtained a global connected set of solutions for a general class of vorticities. 
In the infinite-depth case, the author [Hur06] extended the methods in [CS04] but 
when the vorticity is small, non-negative and monotone with depth. As a matter 
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of fact, the vorticity function of Gerstner's trochoidal waves [Ger09] is not treated. 
The present purpose is to establish for a general class of vorticities a global bifur- 
cation result for Stokes waves on deep water. The main result (Theorem 12. 1[) does 
not have any restriction on the sign of the vorticity. 

Following [CS04], there have been vigorous activities in studies of traveling water 
waves with vorticity. In [KS08b, KS08a], waves found in [CS04] arc numerically 
computed. For arbitrary vorticities, small-amplitude solitary waves over channels 
of finite depth are constructed in [Hur08a] and independently in [GW08]. For 
arbitrary vorticities, symmetry property is studied in [Hur07, CEW07] for periodic 
waves of finite depth and in [Hur08b] for solitary waves. Partial results regarding the 
Stokes conjecture on a "limiting" wave are given in [Var08, Var09]. Hydrodynamic 
stability of periodic waves of finite depth is studied in [HL08]. 

The main result Theorem 12.11 states that for a general class of vorticities a global 
connected set of nontrivial Stokes-wave solutions exists and that the continuum 
contains a sequence of solutions for which either the speed of wave propagation be- 
comes arbitrarily large (cavitation) or the relative flow speed somewhere in the fluid 
region becomes arbitrarily close to zero (stagnation). Furthermore, Theorem 12.21 
states that if the vorticity is non-positive and monotone with depth then stagnation 
can only occur at the wave crest. If the vorticity is non- negative and monotone with 
depth, Theorem 12.31 states that in the event of stagnation, it occurs either at the 
infinite bottom or on the free surface. For non-negative and monotone vorticities, 
if the relative flow speed is bounded along the continuum and if vorticity is, in 
addition, sufficiently small then stagnation occurs either at the infinite bottom or 
at the wave crest. 

The class of vorticity functions admissible in Theorem 12.11 extending consid- 
erably that in the earlier work [Hur06], does include the vorticity distribution of 
Gerstner's trochoidal waves [Ger09]. Furthermore, for small, non-negative and 
monotone vorticities, Thcorcm l2.3l improves the result in [Hur06] in that the small- 
ness condition of the vorticity (|2.7p is more straightforward than that in [Hur06] . 

Cavitation is considered to be unphysical. In fact, in the irrotational setting 
[Tol96,BT03] as well as in the finite-depth case with vorticity [CS04], speed of wave 
propagation is a priori bounded and hence a "limiting" wave along the continuum 
must exhibit stagnation. An important open problem is to obtain bounds for speed 
of wave propagation in the infinite-depth case and with vorticity. In the irrotational 
setting, a bound for traveling speed [Tol96, BT03] is established by studying the 
kernel associated to Nckrasov's integral equation. Unfortunately, such an integral 
representation of solutions [Nek51,Tol96] critically hinges upon irrotationality, and 
presently no analog in the rotational setting is available. Section [6] presents the 
reformulation of the problem with vorticity, which is potentially useful in obtaining 
an integral representation of solutions. 

Another important open problem is to show for non-positive vorticities that 
stagnation occurs at the wave crest, i.e., to remove the monotonicity assumption 
from Theorem 12. 21 A difficulty lies in that the relative flow speed does not possess 
a maximum principle. In the finite-depth case, the result is established in [CS07, 
Var08] by studying the maximum of the relative horizontal speed along the free 
surface, along the bottom and below the wave crest. 
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Ideas of the proof The present treatment is influenced by [CS04, Hur06], based 
on the reformulation via a partial hodograph transform of the original semilinear 
elliptic boundary value problem with free surface (Section 12. ip as a quasilincar 
elliptic boundary value problem in a fixed semi- infinite strip (Section 12. 3p . In 
the finite-depth case [CS04], the Fredholm property of the operator describing the 
problem follows in the standard manner by the Schauder theory for elliptic problems 
and the embedding properties of Holder spaces of functions in a bounded domain. 
Consequently, a global connected set of solutions is obtained as an application 
of a generalized degree theory [HS98] and global bifurcation theory [Rab71]. In 
the infinite-depth case, unfortunately, the unboundedness of the domain prevents 
the operator from being Fredholm, and thus a degree-theoretic argument is not 
directly applicable; see Section l2~4l It is noteworthy that in the irrotational setting 
(of infinite-depth) [Tol96, BT03], the problem further reduces to an equation for 
a quantity defined at the one-dimensional free surface, and specialized theory of 
bifurcation applies. 

In order to overcome the failure of Fredholm property, as is done in [Hur06] , the 
operator is approximated by a sequence of Fredholm operators. The framework 
of a generalized degree theory and global bifurcation theory then applies to each 
approximate problem, and a global connected set of its nontrivial solutions is con- 
structed; see Section [31 Section |4~T1 and Section |4~21 The heart of the matter of the 
proof is in Section 14.31 to take the limit of the continua of approximate solutions 
and to show that the limit set of nontrivial solutions of the original problem is con- 
nected and unbounded in the Banach space in use. The connectedness of the limit 
set necessitates, in light of Wyburn's lemma ( [Why58], see also Theorem I4.5[) . a 
uniform decay at infinity of nontrivial solutions of the original problem. In the ear- 
lier work [Hur06] , ad hoc arguments enforce that for a restricted class of vorticitics 
solutions decay exponentially at infinity. Here, an exponential decay of solutions is 
established in Lemma 14.71 by consideration of a Phragmcn-Lindclof theorem, and 
it does not require any special property of the vorticity. Then, in Lemma 14.101 ro- 
bust nodal properties along the continuum assert that the continuum of nontrivial 
solutions of the original problem is unbounded. 

2. Formulation and the main result 

A detailed account is given of the passage from the traveling-wave problem of 
(|1.1|) - (|1.3|) to an abstract operator equation in a Banach space. The main results 
are stated. The failure of the Fredholm property of the operator is discussed, and 
approximate problems are designed. 

2.1. The vorticity-stream function formulation. The traveling- wave problem 
of (jl.ip ~ p.3p seeks for a solution, for which the wave profile, the velocity field and 
the pressure have the space-time dependence (x — ct, y), where c > is the speed 
of wave propagation. In the frame of reference moving with the speed c, the wave 
profile and the flow underneath it appear to be stationary. Let 

Cl v = {(x, y) : — oo < x < oo, — oo < y < v( x )}> S r] = {(x,rj(x)) : — oo < x < oo} 

denote, respectively, the (stationary) fluid domain and the free surface. 

In studies of traveling water waves, it is customary to introduce the (relative) 
stream function tjj(x,y), defined fl v as 

(2.1) ip x = -v, i' y =u — c 
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and -0(0, 77(0)) = 0. Accordingly, we formulate the traveling-wave problem of (jTTTJ)- 
(|1.3|) as the free boundary problem as: 

for a function 7(7") defined for r € [0, 00) and for a parameter c > 0, find a curve 
y = i~i(x) defined for x <E K and a function ijj(x,y) defined in Q v such that 



(2.2a) tp v <0 in fi^, 
and 

(2.2b) -Atp =-y(V») in f2„ 

(2.2c) V=0 on S„ 

(2. 2d) |VV>| 2 + 2 5 y = on S v , 

(2.2e) VV' — > (0, — c) as y — > —00 uniformly for x. 



The derivation of (|2.2[) is detailed in [Hur06, Section 2]. 

The condition (|2.2a|) means that no stagnation poini[| exists in the fluid region. 
Field observations [Lig78] as well as laboratory experiments [TK97] indicate that 
for wave patterns which are not near the spilling or breaking state, the speed of 
wave propagation is in general considerably larger than the horizontal velocity of 
any water particle. 

The no-stagnation condition (|2.2a|) guarantees that [Hur06] the vorticity is glob- 
ally a function of the stream function, denoted by lu — "f(ip). It is reasonable to 
require that 7(7") — > as r — > 00. Furthermore, the function 

(2.3) r(p)= [%(- P 'W 

Jo 

is required to be bounded for —00 < p ^ 0. Let 

p—00 

r inf = inf F(p), Too = / 7(-p)dp. 

-oo<p^0 J Q 

If the vorticity is non-negative and monotone with depth, i.e. if j(r) ^ and 
7 ; (r) ^5 for r e [0,oo), then (|2.2a|) is redundant. Indeed, by the maximum 
principle and the Hopf boundary lemma, any solution -0 of (|2.2b ) subject to (12. 2b ) 
must acquire (|2.2a|) . 

The boundary condition (|2.2b ) means that the free surface itself makes a stream- 
line, while (|2.2d ) is a manifestation of Bernoulli's law which states that the quantity 
I 1 2 + 2 gy is a constant on the free surface. The Bernoulli's constant only serves to 
relocate the origin in the y-direction and by adding an arbitrary constant to (|2.2H ) 
changes neither the free surface nor the velocity distribution in the fluid region. 
Thus, without loss of generality, the constant is taken to be zero. The hydrostatic 
pressure in the fluid region is given by 

(2.4) P(x, y) = P atm - i| y)\ 2 -gy + T(-^{x, y)). 

In view to the Stokes wave problem, (|2.2p is further supplemented with the 
periodicity and symmetry conditions that rj(x) and ip(x, y) are even and 2L-periodic 
in the x-variable, where 2L > is the wavelength. 



* By a stagnation point we mean a point where ip y = 0. It is a slight abuse of terminology, 
since traditionally |V»/i| = at a stagnation point. 
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In this setting, L and c are considered as parameters whose values form part 
of the solution. The wavelength L, in existence theory, is independent of other 
parameters, and hence it is held fixed in the sequel. On the other hand, the speed 
of wave propagation c serves as the bifurcation parameter and it varies along a 
solution continuum. 

In case 7 = 0, namely in the irrotational setting, (|2.2b ) reduces to the Laplace 
equation and the nonlinearity of the problem resides only at the free surface. A 
nontrivial vorticity, in stark contrast, introduces additional nonlinearity in the field 
equation (|2.2b ). and it significantly complicates analysis. 

2.2. The main results. For a nonnegative integer k and for a € (0, 1), a domain 
ft in R 2 is called a C k+a domain if each point on its boundary, denoted by dfl, 
has a neighborhood in which dVl is the graph of a C k+a function. Given a C k+a 
domain fl in the (x, y)-plane (not necessarily bounded), we define 

(2.5) C k + a (Tl) = {/ e C k+a (ty : f is even and 2i-periodic in the x-variable}, 
where C fc+Q (f2) is a Holder space under the norm 

k 

ll/llc*+«(n) = E ™P\d k y 2 d k x if(x,y)\ 

k 1 +k 2 =0 n 

\d k ^d k -f{x,y)-d k -d k -f{x',v')\ 

fc+SU ((x - x') 2 + (y - y') 2 ) a/2 ' 

n 

This notation is extended in an obvious way to the case when a = and to functions 
of a single variable. 

The main result of this article concerns the existence of nontrivial Stokes waves 
on deep water for a general class of vorticities. 

Theorem 2.1. Let L > be held fixed. Suppose that the vorticity function 7 € 
C 1+Q ([0,oo)) ; a € (0, 1), satisfies that 7(7-) G 0(r~ 2 ~ 2p ) as r — > 00 for some p > 
and that 

(2.6) j U{2T(p) - 2r inf f 2 + (£) 2 (2r(p) - 2r inf ) 1/2 ) e 2 ?dp < g. 

There exists a connected set C of solution triples (c, 77, ip) in the space K + x 
C*3+ Q (R) x C$+ a (Ti v ) of the system (j2^bV ((2~2b ) such that 

(i) C contains a trivial solution which corresponds to a horizontal shear flow 
ip x = under the flat surface 77 = and 

(ii) there is a sequence of solution triples {(cj,r)j,ijjj)} inC, for which 

either lim Cj =00 or lim sup d y ipj(x,y) = 0. 

Moreover, each nontrivial solution triple (c, 77, ip) in C enjoys the following proper- 
ties: 

(iii) i] and "ip are even and 2L-periodic in the x-variable; 

(iv) 1] has a single maximum (crest) at x — and a single minimum (trough) 
at x = ±L per wavelength; 

(v) the wave profile is monotone from crest to trough, i.e., 7] x {x) > for —L < 
x < and rj x (x) < for < x < L; 
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(vi) the speed of wave propagation is larger than the horizontal particle velocity 
everywhere in the fluid region, i.e. ip y < in the fluid region; and 

(vii) Tp x (x,y) > for —L < x < and ip x (x,y) < for < x < L. 

The condition (|2.6|) ensures local bifurcation. A more general condition for local 
bifurcation is in (|4.6| . If T is small, then (|2.6|) is valid. 

Theorem 12.11 presents two alternatives in (ii). If the first alternative realizes, 
the speed of wave propagation along the continuum becomes unboundedly large, 
and correspondingly, the hydrostatic pressure becomes unboundedly low. This 
phenomenon is called cavitation. The second alternative means that the continuum 
contains waves whose relative flow speed somewhere in the fluid region becomes 
arbitrarily close to zero. In other words, there is a region of almost stagnant fluid, 
a region carried along by the traveling wave. This phenomenon is called stagnation. 

If the vorticity is monotone with depth, then the conclusion in (ii) can be refined. 

Theorem 2.2. Under the hypotheses of Theorem \2.1\ if, in addition, 7(7") ^ 
and j'(r) ^ for ^ r < 00, then the conclusions of Theorem ] 2. 1\ holds with (ii) 
replaced by 

(ii') there is a seguence of solution triples {(cj,T)j,ipj)} C C, for which 
either lim Cj —00 or lim dyipj(0,r]jO)) = 0. 

j — >oo j — >oo 

Theorem 2.3. Under the hypotheses of Theorem \2.1[ if, in addition, 7(7") ^ 
and j'(r) ^ for ^ r < 00, then the conclusions of Theorem ] 2. 1\ holds with (ii) 
replaced by 

(ii") there is a seguence of solution triples {{cj,r]j,ifjj)} C C, for which 
either lim Cj —00 or lim Cj = or lim max d y ipj(x,r]j(x)) = 0. 

j — >oo j — >oc j — >oo O^x^L 

If ^-> y (±L,ri(±L)) < M for all (c,n,ip) G C. and if, in addition, 7(0) is suffi- 
ciently small so that 

(2.7) 3 + 7 (0)AO0, 

then (ii") is further replaced by 
(ii'") there is a seguence of solution triples {(cj,r]j,ipj)} C C, for which 

either lim Cj =0 or lim d y ijjj(0,r)j(0)) = 0. 

j->oo j— >oo 

Theorem 12 . 21 states that if the vorticity is non-positive and monotone with depth 
then stagnation, if occurs, must be at the wave crest. If the vorticity is non-negative 
and monotone with depth, Theorem 12.31 states that when the second alternative in 
(ii) of Theorem 12.11 realizes, stagnation occurs either at the infinite bottom or 
somewhere on the free surface. If the relative flow speed at the wave trough is 
bounded along the continuum, and if, in addition, the vorticity is sufficiently small, 
then cavitation does not occur and stagnant occurs either at the infinite bottorm 
or at the wave crest (see (ii'") in Theorem 12. 3p . 

The conclusion of Theorem 12.11 or Theorem 12.21 in case of zero vorticity, partly 
recovers the well-known result ( [Tol96] , for instance) that the continuum of irrota- 
tional Stokes waves contains a "limiting" wave with stagnation at the wave crest. 
In other words, the second alternative in (ii) of Theorem 12 . 1 1 or (ii') of Theorem 12. 21 
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occurs. Theorem 12.11 or Theorem 12 . 2 1 also explains the existence of Gerstner's tro- 
choidal waves [Ger09]. Indeed, the vorticity function corresponding to Gerstner's 
waves [OS01], given by 

e 26(^) 

where ^ m < 1 and — oo < b < 0, is non-positive and monotone with depth. 

Theorem 12 . 31 improves the result in [Hur06]. Cavitation, if occurs, is shown to be 
a consequence of that the speed of wave propagation becomes unboundedly large. 
Furthermore, (|2.7[) gives a more straightforward smallncss condition of the vorticity 
than that in [Hur06]. 

In the finite-depth case [CS04], the second alternative in (ii) of Theorem 12.11 
(stagnation) realizes. Instead of speed of wave propagation, in [CS04] Bernoulli's 
constant serves as the bifurcation parameter, and the parameter values are shown 
to be subcritical; solitary water waves of small-amplitude bifurcate for supercritical 
values of parameter [Hur08a] . While solitary water waves are not expected to exist 
in the infinite-depth case [Cra02,Hur09], nevertheless, I conjecture that a limiting 
Stokes wave with vorticity exhibits stagnation. 

In the irrotational setting [Tol96, BT03] speed of wave propagation of Stokes 
waves on deep water is shown to be a priori bounded by studying the kernel asso- 
ciated to Nekrasov's integral equation. It is noteworthy that even when existence 
theory is based on Babenko's pseudo-differential equation [BT03], a bound for trav- 
eling speed uses Nekrasov's integral equation. Such an integral representation of 
solutions, unfortunately, critically depends on that the stream function is harmonic, 
and thus it is not readily available in the rotational setting. Section [5] presents the 
reformulation of the problem via a quasi-conformal transform, which has structural 
similarity to the formulation of the irrotational problem in [LC25], and thus it is 
potentially useful to obtaining an integral representation of solutions. 

In the finite-depth case [Var08, Var09], if the vorticitiy is non-positive (not neces- 
sarily monotone with depth), stagnation is shown to occur at the wave crest. The 
same result is expected to hold in the infinite-depth case, but no proof is given 
presently. Section [5] collects properties of Stokes waves of infinite depth which are 
relevant to study the location of stagnation points. 

2.3. Reformulation: reduction to an operator equation. Under the no- 
stagnation condition (|2.2ap , exchanging the roles of the y-coordinate and ifj offers 
a reformulation of (|2.2b )- (|2.2b ) in a fixed domain, which serves as the basis of the 
existence theory. 
Let 

q = x and p = —ip(x, y) 

be new independent variables. They map the fluid region of one period {(x, y) £ 
: — L < x < L} to the fixed semi-infinite strip {—L,L) x (— oo, 0) in the {q,p)~ 
plane and the free surface of one period {(x,rj(x)) : — L < x < L} to the top 
boundary (— L,L) x {0} of the strip. Let 

R = {{q,p) : -L < q < L, -oo <p < 0}, T = {(q,0) :-L<q<L}. 

Accordingly, the depth function 

h(q,p) = V 
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(2.8) h q = -02., h p 



replaces the dependent variable. It is straightforward to show that 

h = \_ 

By the above partial hodograph transform, the semilinear elliptic free boundary 
problem (|2.2[) is reformulated as the following quasilinear elliptic boundary value 
problem in the fixed domain R: 

dUJi) (1 + h 2 q )h pp - 2h p h q h pq + h 2 p h qq = -j(-p)hl in R, 
flHfc >) 1 + 2ghh 2 p + h 2 q = on T, 

(|2.9b ) V/i = (h q , h p ) — > (0, 1/c) as p — > — oo uniformly for g, 



where /i p > in i? and /i is even and 2L-periodic in the g-variable. 

It is established in [Hur06, Lemma 3.1] that the above formulation is equivalent 
to f2~2]) . 

A preliminary step of obtaining an operator equation for (|2.9p is to identify its 
trivial solutions. 

Lemma 2.4 (Trivial flows). Given 7 <G C 1+Q ([0, 00)), a G (0,1), for each A G 
(— 2rj n f,oo) i/ie system <\2.9\i has a solution 



(2.10) MP) 



-(p;A)= [ P (\ + 2r(p'))-^dp>-±, 
Jo 2 9 



which corresponds to the shear flow in the horizontal direction 

(My(y) = -(x + 2r( P (y))y/ 2 

under the flat surface r]{x) = 0, where p(y) is an inverse of (|2.10p . 
The proof is in [Hur06, Lemma 3.2], and it is omitted. 

In the bifurcation analysis below, instead of c the square of the (relative) up- 
stream speed A = (h' tr )~ 2 (0) = (^t r ) 2 (0) of the trivial flow (|2.10[) serves as the 
bifurcation parameter. For each A G (— 2Fi n f , 00), the speed of wave propagation is 
determined by A by c 2 = A + 21^00 . 

It is convenient to make use of the shorthand 

(2.11) a(A) = a( J9 ;A) = (A + 2r(p)) 1 / 2 . 
The derivatives of h tr can be expressed in terms of a as 

h' tr (p) = a-^p-X), h'Up) = - 7 (™P)«- 3 (P;A). 
Note that a(A) is bounded for each A G (— 2Fi n f , 00). 

In order to tackle the existence question for solutions of (|2.9p via bifurcation 
theory, we need to further reformulate the problem as an abstract operator equation 
in the form F(X,w) = 0, where w belongs to a Banach space. To this end, let 

(2.12) h(q,p) = h tr (p)+w(q,p). 

Then, w(q,p) — > as p — > —00 uniformly for q. 
We introduce the function spaces in use. Let 

X = { f G C^+ a (R) ■ dfrd^w G o(l) as p -» -00, ki +k 2 ^3 uniformly for q }, 

Y x = { f G Cl+ a {R) ■ Bfcefrw G o(l) as p -» -00, h + k 2 < 1 uniformly for q}, 
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and Y 2 = C 2 + a {T). Recall that the subscript per means evenness and 2L-periodicity 
in the g-variable. Let Y = Y\ X Y 2 with the product topology. We equip X and Y 
with the Holder norms (thus rendering them Banach spaces): 

|| ' \\X ■= || ' Ilc3 + °("R) ' II ' 11^ := II ' 11^1 + II ' WY2 : 

where || • \\ Yl = || • and || • ||y 2 = || • ||c 2 +°(t)- Let Z = C° er (R) have the 

usual maximum norm || • \\z = \\ ■ \\c°(r)- 

The operator form of the Stokes wave problem is then given 7(7") defined for 
r £ [0,oo) to find a nontrivial solution (A,u>) € K x X of 

(2.13) F(X,w)=0, 
where 

(2.14) F(X, w) = {F 1 (\,w),F 2 (X,w)) : R x X -> Y, 



(2.15) 



Fi(X,w) =(1 + w q )w pp - 2(a (A) + w p )w q w pq + (a (A) + w p ) w qq 
+ 7(-p)(a" 1 (A) + w p f - 7 (-p)a- 3 (A)(l + w 2 q ), 

(2.16) F 2 (X,w) = [l + (2.gw-A)(A- 1 / 2 +u.p) 



\2 , 2 
W q 



2.4. Approximate problems. In the finite-depth case with vorticity [CS04] as 
well as in the irrotational setting (of infinite depth) [Tol96, BT03], the key to a 
successful existence theory "in the large" for Stokes waves lies in a generalized 
degree theory and global bifurcation theory. 

The rotational Stokes-wavc problem in the finite-depth case in [CS04] takes the 
same operator equation ^(A, w) = as in the infinite-depth case (where F is in 
(|2.15[) and (|2.16[) ) but with the important difference that in the finite-depth case w 
is considered in a finite rectangle, whereas in the infinite-depth case it is considered 
in the semi-infinite strip R. 

In the finite-depth case, the equation F(X, w) = gives an elliptic boundary value 
problem in the bounded domain, and the Fredholm property of F follows [CS04] 
from the Schaudcr theory for elliptic problems and the compact embeddings of 
Holder spaces of functions in the bounded domain. The existence of Stokes wave 
solutions in the finite-depth case [CS04] then uses the generalized degree theory, 
adapted by Healey and Simpson [HS98] for a general class of nonlinear Fredholm 
operators, and global bifurcation theory [Rab71]. 

In the infinite-depth case, unfortunately, a similar approach fails. Denoted by 
F W (X, w) is the Frechet derivative of F in the second argument at (A, w) € M x X. 
A straightforward calculation yields that 

F w (X,w) = (A(X,w),B{X,w)), 

where 
(2.17) 

A(A, w)[(p] =(1 + w 2 q )ip pp - 2(aT 1 (A) + w p )w q ip pq + (a _1 (A) + w p ) 2 ip qq 

+ (-2w q w pq + 2(aT 1 (A) + w p )w qq + 37(-p)(aT 1 (A) + w p ) 2 ) ip p 
+ (2w q w P p - 2(a~ 1 (A) + w p )w pq - 2"f(-p)a~ 3 (X)w q ) ip q , 

(2.18) 

B(X, w)[<p] = \2{2gw - A)(A" 1/2 + w p )tp p + 2w q Lp q + 2g(\~ 1/2 + w p ) 2 (p . 
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We shall show in Lemma l3~Tl that the closed-ness of the range of F W (X, w) : X — ► Y 
is equivalent to the unique solvability of its "limiting" problem 

<Ppp + (A + 2T oc )~ 1 ip qq = 

in the infinite strip {(q,p) ■ —L < q < L, — oo < p < 00} in the Cp er class; see 
also [VV03]. However, the spectrum of the operator d 2 + (X + 2T oc )~ 1 d 2 defined in 
the infinite strip is (—00, 0], and the limiting problem has infinitely many solutions. 
In the infinite-depth case the operator F defining the Stokes-wave problem is not 
Frcdholm, and (generalized) degree theory may not be directly applicable. 

This difficulty can be overcome by studying a sequence of "approximate" prob- 
lems 

(2.19) F e (X,w)=0, 
where e > and 

(2.20) F e (\, w) := (F 1 (X,w)-ew,F 2 (X,w)). 

Then, the Fredholm property of F e for each e > follows from that its limiting 
problem 

<P P p + (A - 2T oo y 1 ip qq - eip = 

in the finite strip {(q,p) ■ —L < q < L, —00 < p < 00} admits only the trivial 
solution; see Lemma 13. II and Lemma 13.31 



3. Generalized degree for the approximating operators 

For S > let us define the set 
(3.1) 

{2A S 
(A, w) G R x X : X > -2T inl + 5, a _1 (A) + w p > S in R, w < — — on T 

The purpose of this section is for each 6 > and for each e > to establish several 
properties of the operator F e on the set Os needed to define a topological degree. 

First, for each 5 > and e > by continuity of F £ it follows that Os is open in 
1x1. It is straightforward that 

F^(X,w) := (A^X 7 w),B(X,wj):X^Y, 

where ^4 c (A,u>) = A(X,w) — ew and A(X,w) and B(X,w) are given in (|2.17p and 
(|5.13p . is continuous. Furthermore, F e : R x X — > Y is at least twice continuously 
Frcchet differcntiable. 

Next, the principal parts of operators A t (X,w) and B(X,w) are denoted by 

(3.2) A P (X, w)[ip] = (1 + w q )ip pp - 2(a~ 1 (A) + iv p )w q ip pq + (a _1 (A) + w p ) 2 ip qq , 

(3.3) B p (X,w)[ip] = 2(2gw - X){X~ 1/2 + w p )<p p + 2w q ip q 

respectively. For each 8 > and for each (A, w) G Os note that the differential op- 
erator A e (X, w) is uniformly elliptic with coefficient functions bounded in C 2+a (R); 
the coefficients of the principal part satisfy 

4(1 + ^)(a- x (A) + w p ) 2 - A{a-\X) + w p ) 2 w 2 > AS 2 . 
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Also, note that the boundary operator B(X,w) is uniformly oblique in the sense 
that it is bounded away from being tangential; the coefficient of ip p in B(X,w) 
satisfies 

\2(2gw- A)(A~ 1/2 +w p )\ > S 2 on T. 

For each 8 > and for each (A, to) e Os thus F^(X,w) = (A e (A, to), B(X, w)) 
satisfies the following Schauder estimate [ADN59] 

(3.4) \\tp\\ x < C(\\A*(X,w)[tp}\\ Yl + \\B(X,wM\\ Y2 + |MU) 

for all tp E X, where C > is independent of tp. 

Recorded in the next lemma is the Fredholm property of F^(X,w). 

Lemma 3.1 (Fredholm property). For each 8 > and for each e > 0, for each 
(A, w) 6 Os the linear operator F^(X,w) = (A e (X,w), B(X,w)) : X — > Y is a 
Fredholm operator of index zero. 

Proof. The first step is to show that F£,(\,w) : X — > Y is semi-Fredholm. That is, 
its range is closed in Y and its kernel is finite-dimensional. 

Let {(fj} be a bounded sequence in X. Let a sequence {(j/y,J/2j)} converge to 
(j/i, j/2) in y as j — ► 00, and let 

F^(X,w)[(pj] = {A e (X,w),B{X,w))[p 1 } = (yy.ifty), j = 1,2,.... 

It is immediately that y?j — * tp in Cp er {R') for some </? for any bounded subset R' 
of R. By continuity, then 

F wi\w)[p] = (j/1,2/2). 

Our goal is to show that tpj — > 93 in X. 

We claim that y>j — ► <p in Cp er (R). Suppose, on the contrary, that there exists a 
sequence {(<Zj,Pj)} in R such that pj — > —00 as j — > 00, yet 

(3.5) IVifehft') - ^fe,Pj)l ^ k > for all j 
for some ft. For each j, let us form the "shifted difference" 

&j(q,p) = <Pj[Q,p+Pj) - <p{<i,p+Pj) 

defined in the shifted domain Rj := {(q,p) € R : —00 < p + pj < 0}. By construc- 
tion, dj satisfies 

A e j{X,w(;- +Pj))[tij] = Vij(-,-+Pj) - Vi(-r+Pj) in Rj, 

where the operator Aj(X, to) is obtained by shifting the coefficient functions of 
A e (X, to) by — in the p-direction, i.e., by replacing a(p; X) by a(p + pj]X) and 
7(— p) by 7(— p — pj); the value at (g,p) of the function to(-,- +Pj) is given by 
w ( a iP + Pi), Uij{') ' ~^Pj); J/i('i ' + are defined in the same manner. 

Passing to the limit as j — > 00 of the above, we obtain that the (pointwise) 
"limiting" function $0 of $j is in the C° class is defined in the "limiting" domain 

Ro = {(q,p) ■ —L < q < L, —00 < p < 00} 

of Rj and that it satisfies the "limiting" equation 

(3.6) (tfo) PP + (A + 2r oo r 1 (tf ) w -«? = in.Ro. 

The limiting equation is obtained by taking the (pointwise) limit of the coefficient 
functions of Aj(X,w(-,- +Pj)) and yij(-,- +Pj), 2/i(v as j — > 00 and it uses 

that \7w(q,p + pj),\7 2 w(q,p + Pj) — > as j — > 00 for all (</,p) € R and that 
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a(X,p + Pj) = (A + 2FOO) 1 / 2 and j(— p — pj) — ► as j — > oo for all — oo < p < 0. 
Moreover, since 99 is even and 2L-periodic in the q- variable, so is "$o- 

It is standard that (|3.6p admits only the trivial solution do — 0. Indeed, multi- 
plying the equation by $0 an d integrating over i?o yields that 



(3.7) / / ((tfo)p + (A + 2r oo )- 1 (7? )^ + e(^o) 2 ) dqdp = 0. 

This, however, contradicts (|3.5[) . and thus proves the claim. 

Since A e (A, w) is uniformly elliptic with coefficient functions bounded in C 2+a (R) 
and since B(X,w) is uniformly oblique, an application of the Schauder estimate 
[ADN59] yields that 

|| - <p\\x < C(||2/ij - 2/i || Vi + Nstej - 2/2 1 1 v 2 + ||Pj - <p\\z) 

for all j, where C > is independent of tfj and ^. By the above claim, the last 
term of the right side vanishes as j ' ^ 00 . Since the first two terms of the right side 
decreases to zero as j — > 00 by hypothesis, it follows that ipj — > <p in C 3+Q (i?) as 
j ^ 00. That means, the range of F^(X, w) is closed in Y. 

Repeating the above argument for {yij,y2j) = (0,0) then yields that the kernel 
of F e (X,w) is a finite-dimensional subspace in X. Therefore, F^(X,w) is semi- 
Frcdholm. 

The next step is to show that 

F*(A,0) = (pi + a- 2 (X)d 2 q + '3~/(-p)a- 2 (X)d p - el, [-2X^ 2 d p + 2gX~ 1 l] ) 

is Fredholm of index zero. Let 

U = (d 2 p + (A + 2T OQ )- 1 d 2 - el, [-2X^ 2 d p + 2.9A- 1 /] J 

and let us consider the one-parameter family of operators 

(l-t)L c + tF^(X,Q) :X^Y for te[0,l]. 

Note that L e is obtained by replacing the variable coefficients of F e (X,0) by their 
pointwise limit as p — » —00. It is standard from the elliptic theory (see [Kry96, 
Chapter 3], for instance) that L e : X — ► Y is bijective. In particular, it is a Fredholm 
operator of index zero. Since F^(X,0) is semi- Fredholm from the previous step, it 
follows by the homotopy invariance of Fredholm index [Kat67, Chapter 4] that 
F^(X, 0) is also a Fredholm operator of index zero. 

Finally, since Os is connected, the assertion follows by the continuity of Fredholm 
index [Kat67, Chapter 4]. □ 

For our next preliminary result, we need several notations to describe. The 
domain of the operator A t (X,w) is defined by 

(3.8) D(A e (X,w)) = {(/) e X : B{\,w)[<p] =0}. 

Note that A e (X, w) restricted to D(A e (X,w)) is closed in Yy. The spectrum of 
A e (X, w) is defined by 

(3.9) a{X,w) = {fi G C : A e (X,w) - fil : D(A e (X,w)) -> Yi is not bijective}, 

where A e (X,w), D(A £ (X,w)), and Y\ are complexified in the natural way. If /j, £ 
a(X,w) and ker(A e (A, w) — fil) is nontrivial then [i is called an eigenvalue. An 
eigenvalue \i is said to have finite algebraic multiplicity if 

dimker(A e (A, w) - fil) m = dim ker(A e (A, w) - nl) m+1 < 00 
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for some positive integer m. In this case, dimker(A e (A, w) — pil) m is called the 
algebraic multiplicity of \x. 

Lemma 3.2 (Spectral properties). For each 5 > and e > and for each (A, w) € 
Os with |A| + ^ M, where M > 0, there exists a small constant s > and 

positive constants C\,C 2 such that 
(3.10) 

dIMIx ^ | M r /2 !l(A e (A,«;) - /i/)b]|k + | M r /2 !|B(A,u;)M||y 2 + \n\ {a+1)/2 M Z 

for all ip G X and for all /j, £ C satisfying |arg(/i)| ^ 7r/2 + s and ^ C7 2 ^ 1 
sufficiently large, where a € (0, 1) is the Holder exponent inherent from X and Y . 

Moreover, cr(A, w) possesses only finitely many eigenvalues in the sector \arg(ii)\ ^ 
tt/2 + ,s, each of which has a finite algebraic multiplicity. The boundary operator 
B(X,w) : X — > Y 2 is surjective. 

The proof is in [Hur06, Lemma 4.5]. See also [HS98, Proposition 4.4]. 
Our last preliminary result is the properness of F e . 

Lemma 3.3 (Properness). For each 6 > and e > 0, the nonlinear operator F e is 
(locally) proper on Os, i.e., (F 6 ) -1 ^) n D is compact m R x X for each bounded 
set D in Os and for each compact set K in Y . 

Proof. Let {{Xj,Wj)} be a bounded sequence in D C Os- Let {{yij,y2j)} be a 
convergent sequence in K C Y. Let (yij,y 2 j) — * (2/1,2/2) as j ^ 00 and let 

F e (X j ,w j ) = (yij,y 2 j), j = 1,2,.... 

Our goal is find a subsequence of {(Xj, Wj)} which converges inKxX. 

It is immediate that (possibly after relabling) Xj — > X as j — > 00 in R for some A 
and that uij — ► w as j — ► 00 in Cp er (R') for some iu for any bounded subset i?' of 
i?. Moreover, by continuity, 

^ £ (A, W ) = (j/ 1 ,y 2 ). 

It is convenient to write F e in the operator form as 

(3.11a) F({X,w) = A p (X,w)[w] + fi{X,w) — ew, 

(3.11b) F 2 (X,w) = B p (X,w)[w] +f 2 (X,w). 

Here, A p (X,w) and B p (X,w) are the principal parts of A(X,w) and B(X,w), re- 
spectively; fi{X,w) = j{-p){a' 1 (X) + Wp) 3 - 7(-p)a~ 3 (A)(l + w 2 q ) and f 2 (X,w) = 
X- 1 / 2 (2gw-X)(X- 1 / 2 +w p ). 

The first step is to show that uij — > w in C° er (i2). The proof is very similar to 
that in Lenrma l3.il Suppose the convergence does not take place; there would be 
a sequence {(qj,Pj)} C R such that — > —00 as j — > 00 but for some re, 

(3.12) l^jfeiPj) - w(<7j,Pj)l > « > for all 

For each j, as is done in Lemma 13.11 let us consider the function 

Vj{q,p) = Wjfap+pj) -w(q,p+pj), 
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defined in the domain Rj = {(q,p) ■ —L < q < L, — oo < p < —pj} (same as 
Lemma l3Tj) . It is straightforward to show that each Vj satisfies 

= - (M'iXj-'r,:: -+p 3 )) - Af{X, W (; ■ + Pj ))) W(; -+ Pj ) 

- hj(>^j,Wj{; -+P])) + flj(X,w(-, ■ +Pj)) 
+ ew(; ■ +Pj) +2/lj(v +Pj) - yi(; ■ +Pj) 

in Rj. Here, A P (X,w) is obtained by shifting the coefficient function a(p;X) of 
A P (X, w) by — pj in the p-axis, and likewise, fij(X, w) is obtained by shifting a(p; A) 
and j(—p) by — pj in the p-axis; the value at (q,p) of the function Wj(-, ■ + Pj) is 
given by Wj(q,p + Pj) and u>(-,- +Pj), 2/ij(v +Pj)i an d 2/i(v + Pj) are defined 
similarly. 

Passing to the limit as j — > oo of the above, similarly to the proof of Lemma 1 3. 11 
we conclude that there exist the limiting function vq of Vj in the Cp er class, the 
limiting domain Rq = {{q,p) '■ —L < q < L, — oo < p < oo} of Rj, the limiting 
operator dp + (A + 2T oc )~ 1 d q l — el such that 

Mpp + (A + 2T oo y 1 (v a ) qq - ev a = in R . 

Indeed, Xj — > A as j — > oo and \?w(q,p + Pj) — ► as pj — » — oo for all (<?,£>) G i?. 
Since due to the limiting equation the energy integral (13 . T|) of «o is zero, it follows 
that vo = in i?o- This, however, contradicts (|3.12p and proves the convergence of 
{wj} to w as j — > oo in Cp Pr (R). 

Next, since {u>j} is uniformly bounded under the C 3+Q (i?) norm, an interpola- 
tion inequality (see [GT01, Lemma 6.32] and [Kry96, Theorem 3.2.1]) asserts that 
Wj — > w; in C fc +Q (i?) for any k' + a' < 3 + a. Indeed, for any s > there exists a 
constant C = C(s) > such that 

(3.13) |K - u;|lo-'+°'("R) ^ S IK - Hlc 3 +=fR) + C(s)||xyj- - w\\ c0(R) . 

The final step is to employ the Schauder theory to obtain the convergence of 
{wj} in C 3+a (R). By virtue of the decomposition (|3.1ip . the difference Wj — w 
satisfies 

A p (X j ,w j )[w j -w] = yij - (A p (X j ,w j ) - A p {X,w))[w] 

- (fi(Xj,Wj) - fi(X,w)) +e(wj - w) 
B p (Xj,w j )[w J -w]= y 2j - 2/2 - {B p (Xj, Wj ) - B p (X,w))[w] 

-(f 2 (X 3 ,w 3 )~f 2 (X,w)). 

Since (Xj,Wj) G Os, the Schauder estimates [ADN59] applies to the operator 

( A p (Xj ,w j ) 1 B p (Xj,w :j )) to yield 

(3.14) 

IK -w\\ x ^ C{\\A p (Xj,w 3 )[wj -w]\\ Yl + \\B p (Xj,w 3 )[wj - w]\\y 2 + IK -w\\ z ). 

The result of the first step is that the last term on the right side of the above 
inequality tends zero as j — > oo. Since {wj} C X is bounded, by the interpolation 
inequality (|3.13[) it follows that coefficients of A p (Xj, Wj) are equicontinuous in j. 
Since IKHyi — > ilHIvi as j — » oo, moreover, it follows that 

\\(A p (X j ,w j )-A p (X,w))[w]\\ Yl -0 asj^oo. 
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The convergence 

\\fi(Xj,Wj) - /i(A, -uj) || vi ^0 as j -> oo 

follows by that f± consists of polynomial expressions of w p and w q . These together 
with the convergence of {yij} in Y\ yield that 

|| A p (Xj, Wj)[wj — wJHy! — > as j — * oo. 

On the other hand, the standard Schauder estimates and the embedding prop- 
erties of Holder spaces in the bounded domain T confirm 

\\B p (Xj, Wj)[wj — w]\\y 2 —> as j — > oo. 

By (|3.14j) . therefore, Wj — > w in C 3+a (i?). The assertion then follows since X is a 
closed subspace of C 3+a (R) . □ 

With the properties of F e established above in hand, for each 5 > and for 
each e > we define a generalization of the Leray-Schauder degree due to Healey 
and Simpson [HS98] for F e (X,w), where (A, w) E Os- The detailed development 
is in [HS98, Section 4|. Our interest in degree theory lies in that the degree is 
invariant under homotopy and hence it can be used in global bifurcation theory. 

4. Existence theory for rotational Stokes waves 

Undertaken is the study of global bifurcation for (|2.13|) . For each < e < 1, 
the existence of nontrivial solutions of (|2.19| is established in a neighborhood of 
the trivial solution for a parameter value A £ . Then, for each < e < 1 the local 
curve of solutions extends to a global connected set of solutions of (|2.19p . Finally, 
a global existence theory for (|2. 13[> is obtained via abstract bifurcation theory and 
Wyburn's lemma in topology. 

4.1. Local bifurcation for approximate problems. It is readily seen that for 
each e ^ 0, the points (A, 0), where A £ (— 21^, oo), form the line of trivial 
solutions in R x X. The linearization of F e (X,w) about the trivial solution (A, 0) 
is F^(A,0) = (A e (A,0),B(A,0)), where 

A e (A, 0)M = a- 3 (X)(a 3 (X)^ p ) p + a- 2 (X)^ qq - ep, 

B{X, 0)[<p] = -2A 1/2 V P + 2g\~ 1 <p . 

Since F e is continuously Frcchct diffcrcntiablc and F^(X, 0) is Frcdholm of index 
zero, a necessary condition for bifurcation from a trivial solution (A, 0) is that 
F^(X, 0) : X — > Y is not injective, or equivalently, the boundary value problem of 
the self-adjoint equation 

(4.1a) (a 3 (X)tp p ) p + (a(X)tp q ) q - ta 3 {X)tp = in R, 

(4.1b) A 3/ V P = 9<P on T 

admits a nontrivial solution in X. 

Lemma 4.1 (Bifurcation points). Suppose that 7 S C 1+Q ([0, 00)), a £ (0,1), 
satisfies (|2.6p . 

(i) For each e < 1 there exist a unique X e G (— 2Fi n f, gL/ir — 2Fi n f] and a 
unique (up to constant multiple) nontrivial solution tp e G X to (|4.1D . 

(ii) For each s$ e < 1, X° < A e and A c — * A as e -> 0+. 
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Proof, (i) Let ^ e < 1 be held fixed. In view of evenness and periodicity of <p € X 
in the q- variable, we look for a solution of the form ip e (q,p) = <E> e (p) cos k{ir/L)p, 
where k ^ is an integer. Then, $ € solves the ordinary differential equation 

(a 3 {X)v')' - k 2 {n/L) 2 a{X)v - ea 3 (X)v = 0. 

Let us consider for A € (— 2Li n f,oo) the (singular) Sturm-Liouville problem 



(4.2) 



LH) := -(a 3 (A)v')' + ea 3 (X)v = fi(X)a(X)v for - oo < p < 0, 
A 3 /V(0) = gv(0) 



subject to the boundary conditions v, v — > as p — > — oo. Here, the prime denotes 
differentiation with respect to the p-variablc. Our aim is to find a A e such that 
yu(A e ) = — k 2 (ir/L) 2 is a generalized eigenvalue of (|4.2|) (such that L e v — fia(X e )v 
for some v ^ 0), where k is an integer. 

The proof is based on the variational consideration of (|4.2j) . Let us define the 
Rayleigh quotient 

(4.3) R(X) = R (v; A) = -g — — 

J-oo a ( A ) v d P 

and consider the minimization problem 

(4.4) A e (A) = inf{ R e {v; A) : u G ^((-oo, 0)) and v ^0}, 

where 77 1 ((— oo, 0)) denotes the L 2 -Sobolev space in the p- variable. It is well known 
that L e with the boundary conditions in (|4.2[) has the continuous spectrum [e, oo). 
The Rayleigh principle thus asserts that A e (A) is the lowest (generalized) eigenvalue 
of (|4.2p if and only if A e (A) € (— oo, e). Furthermore, such an eigenvalue A £ (A) is 
simple. Our aim is then to find a A c such that A e (A e ) = —k 2 (n/L) 2 . There may 
be multiple solutions, for instance, corresponding to different values of k. Here, we 
restrict ourselves to finding one for k = 1. Note that A e (A) is a C 1 function of A. 
First, for A e [gL/n + 2r in f,oo) it is straightforward that 

(a(X)(ir/L) 2 v 2 + a 3 (X)(v') 2 + ea 3 (X)v 2 )dp 



2„,2 i „3/ \\f„J\2\ 



^ / (a(X)(n/L) 2 v 2 + a 3 (X)(v') 2 )dp 

J — oo 

/0 M 
a 2 (X)vv'dp > 2g / vv'dp = gv 2 (0) 
-oo J — oo 



for every v € if oo, 0)). The second inequality uses the Schwarz inequality. 
Thus, f? e (A) > -(rr/i) 2 , and in turn A e (A) ^ -(vr/L) 2 . 
Next, provided that (|2.6[) holds, one can show that 



2 



A £ (-2r inf ) ^i? £ (e p ;-2r inf ) 

-.9 + /-oo« 3 (- 2r inf)e 2p ^ + e /^a 3 (-2r inf )e 2 f dp 
/° oo a(-2r inf )e 2 Pdp 

By continuity, then, there exists A e € (— 2Fi n f , gi/77 + 2F; n f] such that A e (A e ) 
-(rr/i) 2 . 
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The uniqueness of A e follows by that A e (A) is a monotonically increasing function 
of A as long as A e < 0. The proof is nearly identical to that in [CS04, Lemma 3.4] 
in the finite-depth case (and when e = 0), and hence it is omitted. 

Next, let $ e G oo,0)) be an cigenfunction of (|4. 2|) corresponding to the 

generalized eigenvalue fi(X e ) = —(ir/L) 2 . Since the eigenvalue is simple, <& e is 
unique up to constant multiple. It follows from regularity theory that $ € is smooth. 
Moreover, since 

(A e + 2r inf ) 1 / 2 ^ a(A £ ) < (A e + 2r sup ) 1 / 2 , 
where T sup = sup_ oo< „ <0 T(jp), the comparison theorem for second-order ordinary 
differential equations [CL55, Chapter 8] asserts that $ e (p) > for all — oo < p < 
and it decays exponentially: 

for some constant K > 0. Therefore, tp e (q,p) = $ e (p) cos(7r/L)q is in X. 

(ii) The assertion follows since A e (A) is continuous and it is nondecreasing func- 
tion of e ^ 0. □ 

It follows as an application of the local bifurcation theorem from a simple eigen- 
value [CR71] that for each < e < 1 there emanates from (A e ,0) a local curve in 
Rxlof solutions to ([2~T9"|) . 

Proposition 4.2 (Local bifurcation for approximate problems). Let L > be held 
fixed and let 7 G C 1+Q ([0, 00)), a G (0, 1), satisfy (|2~^|) . 

For each < e < 1, there exist sq > sufficiently small and a C 1 -curve 

Cf oc = {(Hs),w(s)) GRx X : \s\ < s }, 

where each point (X(s) , w(s)) in the curve is a solution to (|2 . 19|) . 

At s = 0, the solution (A(0), w(0)) — (A e ,0) corresponds to a trivial shear flow 
under the flat surface. At s > 0, the corresponding nontrivial solutions (A(s),w(s)) 
enjoys the following properties: 

A(s) = A e + O(s) as s -> 0, 

(4.5) w(s) = s$ e (p) cos(x/L)q + 0(s 2 ) as s -> 0. 

The proof is almost identical to that in [Hur06, Appendix A], and hence it is 
omitted. 

One may replace the sufficient condition (|2.6[) for bifurcation by a more general 
condition that the system 

f (a 3 (X)v')' - ea 3 (A)i) = - (vr/ L) 2 a(X)v for p G (-00, 0), 
1 A 3 /\>'(0) = gv(0) 



(4.6) 



with v, v' — > as p — * —00, admits a nontrivial solution for some A £ (— 2r; n f , 00). 

Remark (Bifurcation at e = 0). A candidate bifurcation point A and a transversal 
solution ip° of F w (A , )[ip] = 0, when e = in Lemma [4.11 exist for the singular 
problem (|2.13[) . But, even the local bifurcation for (|2.13|) is singular. For, F W (X, 0) : 
X — > Y is not Fredholm, and thus local bifurcation theorem (due to Crandall and 
Rabinowitz [CR71], for instance) is not applicable. Power series methods may yield 
a direct proof of small-amplitude solutions of (|2.13|) , but since the present purpose 
is global existence theory, we proceed with singular theory of bifurcation. 
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For the future reference, recorded here are that A e (— 2L; n f , gL/ir — 2L; n f] is 
the unique solution of 

veffM(-oo,o)) J_ oo a(A)w 2 dp 

and that ip°(q,p) = <5>°(p) cos(n/L)q € X is a unique (up to constant multiple) 
nontrivial solution of F^A , 0)[y] = 0. Moreover, > for all — oo < p < and 

it is smooth and decays exponentially as p — > — oo. Finally, A°(A) is monotonically 
increasing with A as long as A < 0. 

4.2. Global bifurcation for approximate problems. Now the presentation is 
for < e < 1 the existence in the large of solutions to (|2.19[) via global bifurcation 
theory. 

Throughout the subsection, 5 > and < e < 1 are held fixed. Let SJ be the 
closure in R X X of the set of all nontrivial solution pairs (A, w) <G Os to (|2.19l) . 
where F e is given in (j2~2"0)) and Os is defined in [[3TTj) . Let C| C 1 x I be the 
connected component of <SJ containing the bifurcation point (A e ,0) determined in 
Lemma 14.11 The local curve of solutions Cf oc constructed in Proposition 14.21 is 
contained in CJ. 

Recall from Section [3] that for each 5 > and for each e > 0, the generalized 
degree due to Hcaley and Simpson [HS98] is successfully defined to F e in the set 
Os- The following global bifurcation result is then immediate. 

Proposition 4.3. For every S > and for every < e < 1, at least one of the 
following holds: 

(i) C| is unbounded in R x X ; 

(ii) C| contains another trivial solution pair (A,0) with A 7^ A e ; 

(iii) C| meets dOs- 

The proof is almost identical to that of [Rab71, Theorem 1.3] except that here 
the generalized degree defined in [HS98] is used in place of the Leray-Schauder 
degree, and hence it is omitted. 

The remainder of this subsection is devoted to refining the result of Proposi- 
tion 14.31 and to establishing some properties of C| useful in characterizing C| in 
terms of bounded subsets of Os- 

Demonstrated below is how exploitation of symmetry rules out the second alter- 
native from Proposition [4T3J In order to state the result precisely, let us denote the 
nodal set and its boundaries by 

R~ = {(q,p) eR:-L<q<0}, T~ = {(q,0) : -L < q < 0}, 

dRJ = {(-L,0) : -00 <p < 0}, OR' = {(0,p) : -00 < p < 0}. 

First, upon examination of (|4.5| . the following nodal properties 

(4.7) w q >0 in iT LIT", 

(4.8) w qq > on dR[, w qq < on dR~ , 

(4.9) %(-£,0)<0, «> OT (0,0) >0. 

hold along the local bifurcation curve CI- The proof is based on properties of 
$ e cos(7r/i)g in Proposition ^. 2\ and it is detailed in [Hur06, Lemma C.l]. 
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Next, since A e (X, w) is a uniformly elliptic second-order linear partial differen- 
tial operator and since B(X, w) is uniformly oblique boundary operator, the maxi- 
mum principle, the Hopf boundary lemma and the edge-point lemma applies to the 
boundary value problem 

(A e (\,w),B(X,w))[w q ] =0 

in R~ . Consequently, each nontrivial solution w of (|2.19p with (A, w) € CJ \ (A c , 0) 

possesses (|4.7j) . (|4~8j) and 

(4.10) 

either w qq (0, 0) < or w qqp (0 7 0) > 0, either w qq (—L, 0) > or w qqp (—L, 0) < 0, 

unless CJ contains a trivial solution (A, 0) other than (A e , 0). The proof is in [Hur06, 
Lemma C.3]. On the other hand, by the variational characterization of A e , it follows 
that if (A,0) e CJ then A = A e . See [Hur06, Lemma 4.8]. See also Lemma l4~T0l 
Therefore, the nodal properties (|4.7p . (|4.8|) and ()4.10[) preserve along the continuum 
CJ\(A%0). 

Finally, our task is to prove t|4. 9|) assuming that (|4. 10|) holds. By evenness, 
w q (0,p) = w qqq (0,p) = for — oo < p ^ and w qp (0,p) = for — oo < p ^ 
0. Moreover, from (|4.8|) it follows w gg (0,0) ^ 0. Differentiating the boundary 
condition i*2(A,u>) = twice in the q- variable results in 

gw qq (\- 1/2 + w p ) 2 + Agw q (\- 1/2 + w p )w qq + (2gw - \)w 2 pq 

+ (2gw - A)(A _1/2 + w p )w pqq + w 2 q + w q w qqq = 

on T. Evaluated at (0,0), it reduces to 

gw qq (\~ 1/2 + w p ) 2 + {2gw - A)(A~ 1/2 + w p )w pqq + w 2 q = 0. 

Suppose that w qq (0, 0) — so that (|4.10p would dictate w pqq (0, 0) > 0. Then, the 
above equation would reduce to 

(2gw~- A)(A" 1/2 + w p )w pqq = at (0,0). 

Since 2gw — A and A -1 / 2 + w p are nonzero, however, w pqq (0,0) = would hold. 
This proves the assertion by contradiction. 
Let us define the open selH 

(4.11) Af = { w G X : w satisfies (j4~7)) - (|43)) } 

and summarize the results. 

Theorem 4.4 (Global bifurcation for approximate problems). For each 6 > and 
for each < e < 1, the global continuum C| either is unbounded in R x X or 
intersects dOg- Each nontrivial solution lying on C| has the nodal configuration 
gH)-gH) ; i.e., CJ \ (A e , 0) c R x Af. 

The following result describes global bifurcation in terms of bounded open sets 
in 0$, and it is crucial in the next subsection in obtaining a global continnum of 
nontrivial solutions to (|2.13|) . 



§ Note that Il4.7l l- ll4.9l l define an open set in X, while 114.71 1. 1 14. 81 1 and 114.101 1 will not define an 
open set. See [CS07] for a more detailed discussion. 
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Theorem 4.5 (see [AT81, Theorem A6]). For each 5 > 0, let S C Os be a closed 
set with (A, 0) G S and let every bounded subset of S be relatively compact inlx X . 
Let C be the maximal connected subset of S containing (A, 0). Then C either is 
unbounded in K x X or meets dOs if and only if dU fl S ^ 6 for every bounded 
open set U in Os with (A, 0) G U. 

Based on Wyburn's lemma [Why58] in topology, the proof is identical to that 
of [AT81, Theorem A6], except that S is confined in Os C R x X so that C can be 
bounded by intersecting dO$ if it is not unbounded in R x X . 

Corollary 4.6. For each 8 > and for < e < 1 sufficiently small, if U is a 
bounded open set in Os with (A , 0) <E U and U C Os then dU nCj ^ 0. 

The assertion follows once Theorem 14.51 applies to the solution branch CJ of 
(|2.19p . Indeed, the result of Lemma [3.31 says that any bounded subset of 6>J is 
relatively compact in R x X, and the result of Lemma HTlT b) says that (A e , 0) G U 
for sufficiently small e > 0. 

4.3. Global existence of rotational Stokes waves. A global connected set in 
R x X of nontrivial solutions to the singular problem (|2.13p is constructed. 
For each S > let 

S s = { (A,w) g Os :F (A, to) = 0, w ^ 0, w G N 

w q G Oflpr 1 -') as p -» -oo } U {(A , 0)} 

for some p > 0. In other words, Ss consists of nontrivial solutions to (|2.13|) with 
the nodal properties (|4.7|) - (|4.9|) and with the decay condition that w q G 0(\p\^ 1 ^ p ) 
as p — » — oo, plus the bifurcation point (A°,0). Let be the maximal connected 
component of the closure in R x X of Ss containing (A , 0) . Our goal is to show 
that C' s is a continnum of nontrivial solutions of (|2. 13[) with the desired properties. 

In order to apply Theorem l4.5l to Ss and C' s , we need to establish that Ss is closed 
and that every bounded subset of Ss is relatively compact in R x X. The relative 
compactness of Ss requires, as we shall see in the proof of Lemma 14.91 a certain 
uniform control at the infinite bottom of functions in a bounded subset of Ss- The 
following result in the spirit of the Phragmen-Lindelof theorem furnishes a uniform 
decay as p — > — oo of solutions w of (|2 . 13[) bounded under the norm of X . 

Lemma 4.7 (Exponential decay). Let 7 G C 1+Q ([0, 00)), a G (0, 1) satisfy 7(7*) G 
( r -2-2 P ) as r ^ oo. For each S > 0, if (A, to) G 5 5 and if |A| + H|x < M for 
some M > 0, then w q enjoys an exponential decay property 

(4.13) \w q (q,p)\ sS M(2 - e /39 )e CT35 /or all (q,p) G i?, 

where [3 > /or some constant K and the constant a > sufficiently small is 

given in (|4.16[) . 

Proof. Differentiating the equation -Fi(A, w) = in the (7- variable yields that 
=(1 + w 2 q )u) qP p - 2(a~ 1 (A) + w p )w q w qpq + (a _1 (A) + w p ) 2 w qqq 

+ (— 2w q w pq + 37(-p)(a~ 1 (A) + w p ) 2 ) w qp + (2w q w pp - 2~/(-p)a~ 3 (\)w q ) w qq 
= :^4 P (A, ui)[w g ] + bi(Vw, V 2 w)w qp + i) 2 (Vw, V 2 w)w qq 
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in i?, where A p (X,w) is the principal part of A(X,w), and &i(Viu,V 2 u>) and 
fe2(V2,V 2 t«) consist of quadratic polynomial expressions in Vw and V 2 w. Since 
7(— p) G 0(| -p| _2_2p ) and w q € 0(|p|~ 1_p ) as p — > — oo, it follows that 

&i(Vw,V 2 u;),fe 2 (V2, V 2 w) G 0(\p\- 1 - p ) as p ^ -oo. 

Since (X,w) £ Os, the operator A p (X,w) is uniformly elliptic. Moreover, a~ 1 (X) + 
w p > 5 > in i?. It is straightforward that 

\hCVw, V 2 w)|, |6 2 (Vw, V 2 w)| < A'A/ 2 for some K > 0. 

Let v = w q , and we consider the function in the nodal set i? _ = {(<?,£>) G ii : 
— L < q < 0}. It is immediate that w solves in R~ the following elliptic second-order 
partial differential equation 

L[v] := A P (X, w)[v] + bi(Vw, V 2 w)v p + b 2 (p, w pp )vv q = 0. 

By the nodal property that w £ Af, it follows that v > in R~ . By evenness of w, 
furthermore, v(q,p) = for q = or q = —L for all — oo < p ^ 0. 
Let us define the auxiliary function 

(4.14) u{q,p)=A(2-e?*)e**-v(q,p), 

in where positive constants A, and /3, a will be determined in the course of the 
proof. 

It is straightforward that 
(4.15) 

L[u] < Ae ap [a 2 {\ + w 2 )(2-e^ q ) - 2aa( a - 1 (A) + w p )w g e /39 - /^(cT^A) + w P )V 9 

+ a6i(Viy, V 2 w)(2 - e^ 9 ) + /36 2 (Vw, V 2 u;)e' 39 ) - L[«] 
Ae" p (2a 2 (l + M 2 )+2f3aKM 2 + 2oKM 2 - 2 5 2 e Pq + f3KM 2 e l3q ^j . 
Let j3 = max ^1, A 2 ^~ ^ so that 

_/? 2 <SV 9 + p K1 ^i^ = f3 e Pi(-f3S 2 + KM 2 ) ^ -^e-? L KM 2 . 

Subsequently, let us choose a > small enough that 

(4.16) 2o- 2 (l + M 2 ) - 4(3<jKM 2 - -e~ pL KM 2 < 0. 

Then, L[u] < in ii". 

Out task now is to examine it on the boundaries of R~ . At the top boundary 
{(<7, 0) : — L < q < 0}, we have 

u(q,p) = A(2 - e~P q ) - v(q,p) ^ A- M = 

if A = M. Since v(q,p) = on the side boundaries {(<?,p) : 9 = or q = — L, po < 
p < 0}, it follows that 

u(q,p) = M{2-e~ l3q )e ap > for g = or q = -L, -oo < p < 0. 

In summary, L[u] < in the domain R~ and u on the boundaries of 
Since A P (A, w) is uniformly elliptic, and £>i(Vttf, V 2 w), b 2 {Vw, V 2 w) G 0(|p| _1_p ) 
as p — > — oo, by the Phragmen-Lindelof theorem [Gil52], it follows that u ^ in 
R~. That is, 

0<«(<?,p) < M(2 - e^e^ in RT. 
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Repeating the above argument for — v on R + yields an analogous inequality. This 
completes the proof. □ 

The above proof applies to Ff(A, w) = Fi(X,w) — ew = 0, e > 0, mutatis 
mutandis to yields an exponential decay of solutions to (|2 . 19[) , analogous to (14.131) . 

Corollary 4.8. For each 6 > and for each < e < 1 if (A, w) € C\ and if 

\M + IIHI-f < M for some M > 0, then w q enjoys an exponential decay property 

(4.17) \w g (q,p)\ < M(2 - e 0q )e ap for all (q,p) E R 

where (3 and a are as in Lemma \4-.7\ 

The exponential decay of solutions of (|2.13[) in (|4. 13)) establishes the relative 
compactness of S$ in R x X. 

Lemma 4.9 (Relative compactness). Suppose that 7 <G C 1+Q [(0, 00)), a G (0,1), 
and 7(7") £ 0(r~ 2 ~ 2p ) as r — > 00 for some p > 0. For each S > 0, any bounded 
subset of S$ is relatively compact in M x X . 

The nonlinear operator F in Os C M x X is not (locally) proper, as is discussed 
in Section [2.41 Indeed, its limiting problem 

v pp + (\ + 2r oo y 1 v qq = 

in the infinite strip {(q,p) ■ —L < q < L, —00 < p < 00} admits infinitely many 
solutions in the C° class. Here, compactness is established only for the solution set 
of F, not for the operator itself. 

Proof. Let {(\j,Wj)} C Ss be a sequence inlxl with |Aj| + < M for all 

j for some M > 0. Note that, F(Xj,Wj) = for all j. It is immediate that {Aj} 
has a convergent subsequence in E. By possibly relabeling the index, let Aj — * A 
as j — ► 00. Moreover, it is immediate that Wj — > w as j — > 00 for some w in i?' for 
any bounded subset B! of R. By continuity, F(A, io) = 0. Our aim is to show that 
{yjj} has a subsequence converging to w as j —> 00 in X. 

As is indicated in the proofs of Lemma 13.11 and Lemma 13.31 a crucial step in 
showing the convergence of {uij} in X is to obtain the convergence of {wj} in the 
Cp er (R) norm. Since {wj} converges in Cp er (R') for any R' C R bounded, it entails 
to show that {wj} decays as p — * —00 uniformly for j. 

It is convenient to write 

(4.18) Wj(q,p)= f d q w 3 (q',p)dq' + w J {0 1 p). 



Note that {d q Wj} and {^(0,-)} are bounded in C 2 + a (R) and C^+ a ((-oo, 0]), 
respectively. 

Since j(—p) G 0(\p\~ 2 ~ 2p ) and d q Wj € 0(|p| _1_p ) as p — > —00 for some p, and 
since |Aj| + ||wj||x ^ M for all j, the result of Lemma l4~7l states that d q uij decays 
exponentially as p — > —00 uniformly for j. More precisely, 

(g, p) I Ce CTp for all (q, p) £ R 

where C > and cr > depends only on 6 and M. An argument of Ascoli type 
then applies to assert that the first term in (|4.18[) has a subsequence converging in 
Cp er (R). Moreover, by the dominated convergence theorem, 

d q Wj(q'p)dq' — > / d q w(q',p)dq' as j — > 00. 
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Next is to examine the latter term in (|4.18|) . Since d q Wj decays exponentially like 
(|4.13|) as p — > — oo uniformly for j, it follows from the classical gradient estimate 
for elliptic equations [GT01, p. 37] that 

(4.19) \d 2 q w 3 (0,p)\ sc C'e a ' p for aU p e (-oo,0) 

where C > is independent of the index j, and er is the same as in (|4.13l) . When 
restricted on the half-line q = 0, by evenness of wj, the equation Fi(Xj,Wj) = 
reduces to 

%Wj + {a-\X) +d p w j )d^w j + 1 (-p)(a- 1 (X) + d p w 3 ) 3 - 7 (-p) a - 3 (A) = 0. 

Since 7 <G 0(r~ 2 ~ 2p ) as r — * 00 for some p > 0, and <9 2 u>j(0,p) decays exponentially 
like (|4.19|) as p — » —00 uniformly for j, then it follows that 7^(0, p) decays as 
p — > —00 uniformly for j. Again, an argument of Ascoli type asserts that {u>j(0,p)} 
has a subsequence which converge in C°((— 00, 0]). In turn, {wj}, possibly after 
relabeling, converges to w in C°(R). 

The remainder of the proof is nearly identical to that of Lemma l3.31 and we only 
outline the various stages of the proof. 

As is done in Q3.11|) in the proof of Lemma 13.31 we decompose F as 

Fi(X,w) = A p (\,w)[w] + fi(X,w), 

F 2 (X,w) = B p (X,w)[w] + f 2 (X,w), 

where A p (X,w), B p (X,w), f\(X,w) and f2{X,w) are the same as in the proof of 
Lemma 13.31 We may write 

A p (X j ,w j )[w j -w} = -(A p (X j ,w j ) - A p (X,w))[w] - (fi(Xj,Wj) - fi(X,w)), 

B p (X J ,w ] )[w ] - w] = -{B p (X 0lWj ) - B p {X lW ))[w] - {hiX^Wj) - f 2 (X,w)). 

bmce A P (X ,w) is uniformly elliptic and B p (X,w) is uniformly oblique, an esti- 
mate of Schauder type [ADN59] states that 
(4.20) 

\\wj - w\\ x ^ C{\\A p (Xj,Wj)[w.j - w]\\ Yl + \\B p (X J ,w. j )[w j - w]\\y 2 + \\wj - w\\ z ) 

holds, where C > is independent of j. Since {wj} is bounded in X and it converges 
in Cp er (R) an interpolation inequality (see (|3 . 1 3[) in Lemma [3~3l or [Kry96, Theorem 
3.2.1]) asserts that Wj — > w in Cp er (R). Accordingly, 

\\A P (Xj ,Wj)[vjj — ^0 as j — * 00. 

Moreover, by the standard Schauder theory and the embedding properties in the 
bounded domain, 

\\B p (Xj 7 Wj)[wj — w]\\y 2 — ► as j — > 00. 

The Schauder estimate (|4.20[) then dictates that Wj — > w in C 3+a (R). Furthermore, 
by continuity, w € X and F(X,w) = 0. 

Finally, since {d q Wj} decays exponentially as p — > —00, it follows that w q G 
0(|p| _ p ) as p — > —00. By Lemma 14.71 in fact, w g decays exponentially as p — > 
—00. This completes the proof. □ 

Next is to show that Sg is closed. If {{Xj, w)} in S$ converges to (A, w) as j — > 00 
and if w is not identically zero, then, by continuity, (A, w) is a nontrivial solution 
of (|2.13[) . Moreover, by Lemma [4.71 and Lemma [4.91 w q decays exponentially as 
p — > —00. That means, (A,w) S 5,5. Thus, it remains to show that if {(Xj,Wj)} 
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converges to (A, 0) as j — ► oo then A = A . For < e < 1, the analogous property 
ensures that the nodal properties preserve along the global continnum C' s \ (A £ ,0) 
of solutions to ([5715]) , 

Lemma 4.10 (Closedness). For each S > 0, if (A, 0) G Ss then A = A . 

Proof. The proof is similar to that of [Hur06, Lemma 4.8]. Let (A, 0) G 5a and let 
{(Xj, Wj)} C Ss, Wj ^ for each j such that (Xj,wj) — ► (A, 0) as j ' ^ 00 in R x X. 

Let = d q Wj /\\d q Wj\\ C 2+c,my Since every c^w., is not identically zero in R, the 
function Vj is well-defined. It is straightforward that each Vj satisfies 

(A(\ j ,w j ),B(\ j ,w j ))[v j ]=0. 

Since {vj} is bounded in C 2+a (R), it follows that Vj — ► u in C% er (R') for some 1; 
for any bounded subset i?' of i?. Moreover, by continuity, 

(4.21) (A(A,0),B(A,0))H =0. 

We claim that Vj — > u in C 2+Q (i?). The proof is similar to that of Lemma T3. 31 or 
Lemma 1431 and thus we only sketch the outline of the proof. 

First, the result of Lemma T4.7I states that Vj(q,p) decays exponentially as p — > 
—00 uniformly for j. Then, by repeating the argument as in Lemma 14.91 one 
accomplishes that Vj -> jj as j -> 00 in Cp er (R). Since| |uj|| C2 +arR) = ^ or eacn h 
subsequently, by interpolation inequality as in (|3"T3"|) . it follows that Vj — > v as 
j ^ 00 in Cp er (R). Next, it is straightforward that u,- — 1; satisfies 

A(Xj,Wj)[vj - v] = (A(X,0) -A(Xj ,Wj))[v] in R, 
BiX^Wj^Vj - v] = (B(A,0)-B(A J ,w J ))H on T. 

Since A(Xj, Wj) is uniformly elliptic and B(Xj, Wj) is uniformly oblique for each j, 
the Schauder estimates [ADN59] yield 

^ \\ V J - V \\c*+"(R) < <?(lWj>«>i)['*'j - U ]llc«(R) 

+ ll-B^-.WjObj - v]\\ C i+c.(t) + hj -v\\ co(R) ), 

where C > is independent of index j. Since ||%Hc2+<»(H) = f° r eacn j an d since 
A(Xj,Wj) is equicontinuous, it follows that 

||(^(A,0)-A(Aj,u;j))H|| c «,( S) ^0 as j ^ 00. 

Moreover, by the standard elliptic theory in a bounded domain, it follows that 

\\B p (Xj, Wj)[vjj — w]\\y 2 — > as j — > 00. 

Therefore, (|4.22|) proves the claim. Furthermore, u, Vti,V 2 w G o(l) as p -> —00 
uniformly for q and |Mlc?2+a(5) = 1. 

By the periodicity and symmetry consideration, it follows that v = d q <p for some 
ip. Indeed, since each Vj is 2L-periodic in the q- variable and since it is of mean zero 
over one period, that is, 

L 

Vj(q,p)dq = for all p G (— 00, 0), 

by continuity, v is also 2L-periodic in the g-variable and it is of mean zero over one 
period. Thus, the assertion follows, where ip is 2L-pcriodic in the g-variable. 
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Let us write (|4.21|) as 
(4.23) {A(\,0),B(\,0))[d q tp)=0. 

Recalling the notation and the result of the nodal preservation in the proof of 
Theorem l4.41 since Vj > on R~ and Vj = on dR^UdR~ for each j, by continuity, 
dgtp ^ on R~ and d q (p = on dR^ U dR~. Furthermore, since d q ip satisfies the 
second-order elliptic partial differential equation (|4.23|) and since d q (p ^ in R~, 
the maximum principle ensures that d q tp > in R~ . Hence, d q ip may be expanded 
as a sine series 

oo 

d q f(q,p) = ^2<pk(p)sin-k(ir/L)q, 

fc=0 

where cpk — * as p — > — oo for all k. Accordingly, (|4.23|) is written as 

OO 

J2(( a3 Wf'kY - k 2 {n/L) 2 a(X)ip k ) sinkq = in R~ , 

k=0 

oo 

^2(-2VX(p' k (0) + 2g/Xip k (Q)) sinkq = on T. 

k=Q 

In particular, ipi solves the boundary value problem 

(a 3 (AK)' = (n/LfaiX)^ for p 6 (-oo,0), 
A 3/ Vi(0)=m(0), 

subject to boundary conditions , c/?^ — * as p — > — oo. We observe that ipi is 
a solution of the Sturm-Liouville problem (|4.2[) when e = with the generalized 
eigenvalue \i = — (it/L) 2 . 

In view of the definitions p~3]) and ([44]) at e = it follows that A (A) < 
i?°(pi;A) = -(tt/L) 2 . Suppose that A (A) < -(tt/L) 2 ; the minimizer v?» of A (A) 
would be an eigenfunction corresponding to the simple eigenvalue A (A) (such that 
R°((p*; A) = A (A)), and hence ip* would not vanish on p £ (-co, 0). On the other 
hand, 

2 f° 

'flip) = y / d q (p(q,p)sm(ir/L)q dq < for all p £ (-co, 0). 

^ J-L 

This contradicts the orthogonality 




a(X)<f t (p)<f 1 (p)dp = 0. 



Therefore, A (A) = —(ir/L) 2 , and A = A follows from the monotonicity of A° 
(see [CS04, Lemma 3.4]). This completes the proof. □ 

The global existence result for nontrivial solutions of (|2. 13[1 is now immediate 
and it is described in the next theorem. 

Theorem 4.11 (Global bifurcation for (|2.13p ). For each S > let C' s denote the 
maximal connected component of the closure m R x X of Ss containing (A°,0). 

(i) The continuum C' s either is unbounded in R x X or intersects dOs ■ 

(ii) Each nontrivial solution lying on C' s has the nodal properties (|4.7|) - (|4.9|) . 
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Proof, (i) By virtue of Theorem l4.5l it suffices to show that if U is a bounded open 
set with (A , 0) G U and U C 0$, then dU fl S$ ^ 0. Let U be such an open set. 

The result of Corollary 14.61 savs that there are sequences {e^} such that €j — ► 
as j — ► oo with < €j < 1 small, and {(A,,-, Wj)} C U such that 

(Xj,Wj) G dU fl for each j. 

It is immediate that (Xj ,Wj) G Os forms a bounded sequence in K x X and that 

(4.24) (Fi(\j,Wj) - ejWj,F 2 (\j,Wj)) = for each j. 

Possibly by relabeling, Xj — > A as j — > oo for some A. It will follow by the methods 
in the proof of Lemma 14.91 that {wj} has a subsequence which converges to w in 
X. By continuity, F(A, w) = 0. 

The remainder of the proof is nearly identical to those of Lemmas 13. 31 or Lemma 
14.91 and we only outline its various steps. 

Without loss of generality, we assume |Aj| + |K||a' < M for each j for some 
M > 0. The result of Corollary 14.81 states that d q Wj(q,p) decays exponentially 
as p — > — oo in (q,p) G R uniformly for j. As is done in Lemma 14.91 then, d^Wj 
decay exponentially as p — ► — oo uniformly for j, and Wj(0, •) decays as p — > — oo 
uniformly for j. In view of (|4.18p . by adapting arguments of Ascoli type, it follows 
that {uij} converges to w as j — > oo in C (R). 

Next, since {wj} is bounded in X, an interpolation inequality similar to that in 
(|3 . 13[) (sec [Kry96, Theorem 3.2.1], for instance) asserts that Wj — > w as j — > oo in 
Cp er (R). Since A p (Xj, Wj) is uniformly elliptic and B p (Xj, Wj) is uniformly oblique 
for each j, Schauder theory [ADN59] asserts that the inequality 

IK - w\\x ^ C{\\A p {Xj,Wj)[wj - w]\\ Yl + \\B p (X J ,w J )[w J - w]\\y 2 + \\u>j - w\\ z ) 

holds for each j, where C > is independent of the index j. It is straightforward 
that 

A p (Xj ,Wj)[wj — w] = tj(wj — w) — (A p (Xj,Wj) — A p (X,w))w 

- (fl(^j, w j) - fl(\w)) + £jW, 

B p (X j ,w j )[w j -w} = -(B p (X ]1 w ] )~B p (X 1 w))w~(h(X ]1 w ] )~h(X,w)). 
By arguments completely analogous to those in the proof of Lemma 13 . 31 then, 
\\A p (Xj,Wj)[wj - w]\\ Yl -^0 as j — > oo, 

and 

\\B p (Xj, Wj)[vjj — w]\\y 2 ^0 as j — > oo 
follow, whence it follows that Wj —> w in X. 

By continuity, F(X,w) =0. If w is not trivial then (A,m;) G C' s , upon recalling 
the definitions of 5,5 and C' s . If w = then by Lemma [4.10i it must follows that 
A = A , and hence (X,w) G C' s . This proves the assertion. 

(ii) The proof, requiring to show the preservation of the nodal property of TV 
along the continuum C' s \ (A , 0), is nearly identical to that of [Hur06, Lemma C.3]. 
Suppose the contrary. Since C' g C S$ is connected, there must be a nontrivial 
solution (A, w) G C' g with w q ^ such that at least one of the nodal properties 
(|4.7|) - (|4.9|) would fail for (A, w). We argue by contradiction using the maximum 
principle, the Hopf boundary lemma, and its sharp form at corner points due to 
Serrin. The detail of the proof is in [Hur06, Appendix C]. This completes the 
proof. □ 
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The proof of (i) in the above theorem entails to take the limit of a bounded set 
in CJ as e — > in the sense to show that a bounded sequence {{Xj,Wj)} for each 
(Aj, Wj) is a solution of the approximate problems (|2.19p with e = 6j converges to a 
solution of the singular problem (|2.13|) as j — > oo if ej — > as j — > oo. The uniform 
decay property of solutions of (|2. 13[) ensures that the convergence takes place in 
the strong topology of K x X . 

The purpose of the following lemma is to obtain bounds for the higher derivatives 
of w in terms of w and w p uniformly along the continuum C' s . 

Lemma 4.12. For each S > 0, if swp^ w j eC < A < oo then swp^,w)eC'(\\ w \\c°(R) + 
\\ w p\\c°(r)) < 00 implies sup (A/!ll)eC , \\w\\ x < oo. 

Proof. Provided that A along C$ is bounded, by the maximum principle it follows 
that w q along C' s is bounded by the maxima in T of w and w p in C' s . Then, by a priori 
estimates of Schauder type due to Lieberman and Trudinger [LT86] for quasilincar 
elliptic partial differential equations with nonlinear oblique boundary conditions, 
it follows that the higher derivatives of w along C' s are bounded by the maximum 
norms in R of w p and w q along C' s . The detail of the proof is in [Hur06, Section 
5.2]. □ 

Remark 4.13. By virtue of Lemma r4.121 in case C' s is unbounded in R x X, there 
is a sequence of solution pairs {(Aj, Wj)} C C' g such that either 

(1) lim Aj = oo or 

(2) Aj is bounded for all j while either lim 1 1 Wj \ \ c0 ^ = oo or lim ||3 p u>j || = oo. 

If the other alternative that C' s intersects dOg realizes, there is a solution pair 
(A, w) £ C' s such that one of the following holds: 

(3) A = 2r inf + (5, 

(4) a -1 (A) + w p = 6 somewhere in R, 

(5) w = 2X 4g S somewhere on T. 

5. Properties of rotational Stokes waves 

Established are properties of solutions of (|2.13|) . and in turn, solutions of (12.21) . 
The main results are proved. 

5.1. Properties of rotational Stokes waves. Throughout this subsection, S > 
is held fixed. For each solution pair (A, if) € C' s of (|2.13p . let (c,r),ip) be the 
corresponding solution triple of (|2.2p . As such, (A, w) and (0,77,1/)) are related via 
the transforms in Section |2~31 More precisely, 

(5.1) c 2 = A-2r oo , r]{x) =w(x,0)- —, 



2.9 



7/1 

(5-2) ^ = - _ ln " , 4> y = — T7TT - T ' 

a L (\) + Wp a L (A)+Wp 

Related to these, 

(5.3) h q =w q = -^-, hp = a~ 1 (X) +w p = — 

Recalled are the notations 

O r) = {(x, y) : —00 < x < 00, —00 < y < i](x)}, S r] = {(x, r/(x)) : —00 < x < 00}. 
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Since (X,w) G K x X it follows that (c,r?, V) el+x C|+ a (M) x C|+ a (n,). 

Established properties are summarized of solutions of (|2.13p . and in turn, solu- 
tions of (|2.2[) . and further properties are inferred. These results are of independent 
interests. 

First, by the nodal properties of w, for any nontrivial solution of (|2.2[) 

(5.4) n x (x) > for - L < x < 0, 

(5.5) ip x (x, y) > for (x, y) G Q,, U SV,, — £ < x < 0, 
and 

(5.6) iM±L, y) > for (±L, y) G fi,, ^(0, y) < for (0, y) G Q„. 
hold. By evenness and periodicity properties of w, furthermore, 

(5.7) r) x (0) = r) x (±L) = 0, 

(5.8) V*(0,y)=V*(±£) = for (0,y), (±L,y) G Q„ 

hold for any solution of (|2.2[) . 

By Lemma T4.71 any bounded solution of (|2.2p enjoys the following exponential 
decay estimate 

(5.9) hM*,lOI <C'e°'y ior (x,y)€fl v , 

where C > large and o 7 > depend only on (5 > 0, c > and the Holder norms 
of rj and ip. By the classical gradient estimate in the elliptic theory [GT01, p. 37] 
it follows that 

(5.10) IV'a:^ (0, j/) | ^ C" e G v as y — > — oo uniformly for x. 

Next, the maximum principle is employed to yield bounds for the relative velocity 
\7ip and the pressure. 

Lemma 5.1 (Bounds for the velocity). A nontrivial solution pair r](x) and tp(x, y) 
of ([221) satisfies 

(5.11) ^MO)) < \V^(x,y)\ 2 -2Y{~^(x lV )) ^ip 2 y (±L,r)(±L)) 
for any (x,y) G tt n . 

The result holds without any restriction on the vorticity function 7. The equality 
holds if the free surface is trivial. 

Proof. The proof is an immediate application of the maximum principle due to 
Spcrb [Spc81, Section 5.2]. See also [Var09, Theorem 3.1]. Here, we include the 
proof for the sake of completeness. 

Let us define a function W : fl v —> M by 

(5.12) W{x,y) = i|W(x,y)| 2 -T(-^(x,y)). 
It is straightforward that 

AW+ ww Wx+ wk Wy=0 ina " 

where L\ = 2 r y(i(})t() x — 2W X and L 2 = 2^(1/))^ — 2W y . Since 

W y = Ipxipxy + Ipylpyy + j{^)4>y = tpxtpxy ~ 4>y4>xx — ► as y -» OO, 
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by the maximum principle and the Hopf boundary lemma, it follows that 
mmW(x,y) < W(x, y) < maxM / (cc,y) (x, y) € £l v - 

Srj S rj 

The assertion then follows since T(tp(x,y)) = on 5, and since \Vip(x,r](x))\ 2 — 
—2gr](x) is nonincrcasing for —L ^ x ^ and is nondccrcasing for x ^ L. □ 

Let us define another function B : fl v — > US by 

(5.13) = i|V#B,y)| 2 +fftf-r(-#i:,y)). 

Note that B(x,y) is the negative of the hydrostatic pressure up to a constant. 

Lemma 5.2 (Pressure estimates). For any solution pair rj(x) and if){x,y) of (|2.2p 
the following inequality 

(5.14) B(x, y) — -max (o, sup j(^))ip ^ m £1, ; 

holds. If, in addition, 7](x) and ijj(x,y) satisfy 

(5.15) g + 7(ip)ipy ^ in £l r; , 
then 

dB 

(5.16) Bix,y) ^ in fL cmrf — — (a;,y) > ok SU, 

on 

where d/dn denotes the outward normal derivative at S v . 

The condition (|5.15|) is valid in the irrotational setting and for non-positive 
vorticities. 



Proof. The proof is similar to that of Lemma 15.11 As in the proof of Lemma 15.11 
it is straightforward that 

under the hypothesis (|5. 15[) . where L\ and L2 are given in the course of the proof 
of Lemma 15. II By the Bernoulli equation (|2.2b ). it follows that B(x,y) = on S v . 
Since 

By = tpxlpxy - Ipylpxx + 9^9 &S V ~* "0°! 

the assertions follow by the maximum principle and the Hopf boundary lemma. 
Repeating the argument in the proof, for an arbitrary vorticity asserts (|5.14|) . □ 

Under the condition (|5 . 1 5[) . since B(x,y) = on the free surface S v it follows 
that 



(l,rjx)-(Bx 7 By) 
while the result of Lemma 15.21 says that 

(- V X,l)-(Bx,By) 



= 0, 

Sr> 



> 0. 
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Thus, B x (x,Tj(x)) < for -L < x < and B x (x,r](x)) > for < x < L. On the 
other hand, for any vorticity function 7, it follows that 



dx \2 Tyy 

= ^ylpxy ~ i>x1pyy = B x (jX, 7]{x)). 

Indeed, ip x + n x ip y = at the free surface S v . Since tp y (x, rj{x)) < 0, it follows that 
ip v (x,rj(x)) is nonincreasing for — L < x < and nondccrcasing for < x < L. We 
summarize this result. 

Lemma 5.3. Under the condition (|5.15|) . any solution pair (r)(x),i/)(x,y)) of (12.2p 
satisfies 

(5.17) ip y (-L,r)(-L))^i> y (x,ri(x))^iP y (0,ri(0))<0 for -L^x^O. 

If the vorticity is non-negative and monotone with depth, then the following 
alternative bound is available for the pressure. 

Lemma 5.4 (Pressure estimate for positive vorticities) . If ~f{r) and Y(r) ^ 
for ^ r < 00 then 

(5.18) B(x, y) + T(-ip(x, y)) < in Sl v . 
Proof. Note that 

B(x,y)+r(-i>(x,y)) = y)\ 2 + gy. 

The (steady) Euler equations yield 

B x = -Ipxtpyy + Tpy^xy and By = —^ylpxx + tpxlpxy + .9, 

whence 

a(b + r(-</>)) + 2 7 '(v.)(s + r(-v>)) = ^4 + i>ly + ^ o in a n . 

Since 7'('0) 0, the maximum principle asserts that B + T(—ijj) attains its maxi- 
mum in £l v either on the free surface or at the infinite bottom. On the other hand, 
since 

d y (B - T(~tJj)) = -ip y ip xx + ipxipxy + .9 - "fWiJy — > g as y — > -co, 

the maximum of B — T(—tjj) is attained at the free surface. The assertion then 
follows since T(—tp) = on the free surface. □ 

As is done for the vorticities which satisfies (|5 . 1 5[) , in the setting of Lemma 15.41 
it follows that 



(l,r) x )-(B x ,B y ) 



= 0, 

Sri 



and 



(5.19) »(*):= (-nx,l)-(Bx,B y )+"f(-r) x ,l)-(Tpx,i> v ) 



> 0. 



Unfortunately, this does not yield the monotonicity of ip y (x,r](x)). Instead, 

■£-tf(x, t)(x)) =2i) y ip xy + 2lj) x {i)xx + 7(0)) 
(5.20) ax , M 

rT ^)+7(0)^(x J ,(x)) 
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If the vorticity is non-positive and monotone with depth, then a simple maximum 
principle yields bounds for ip y , a stronger result than that in Lemma 15. II 

Lemma 5.5. If"f{r) and "f'(r) ^ for r S [0, oo), any nontrivial solution pair 
(rj(x),ip(x,y)) of (|2.2[) satisfies 

(5.21) ip y (0, ?7(0)) ^ ^(x, y) < ^(±L, »y(±£)) 

for any (x,y) E tt v . 

Proof. By differentiating (|2.2b ) in the y-variable, we obtain 

A-0 y = -^f'{ip)ipy in Q v . 

Since ipy < m an d Y( s ) 0j by the maximum principle, ip y cannot have an 
interior maximum in 0^. Moreover, since ip y — > c < as y — > — oo, the maximum 
of ip y must be on S v . The assertion then follows by (|5.17jl . □ 

If the vorticity is non-positive, then the amplitude is bounded by the speed of 
wave propagation. 

Lemma 5.6. If 7(7") ^ and if (rj(x),ip(x,y)) is a solution pair to (|2.2p with the 

parameter c > 0, then 

(5.22) 

< (2 5 ) 3 / 2 (|r,(±£)| 3/2 - \v(0)\ 3/2 ) = \M±L,V(±L))\ 3 l^(0,r;(0))| 3 < cL. 

Proof. By integrating (|2.2b ) in the domain Q~ = {(x,y) € Q, n : —L < x < 0} and 
in the light of the Green's theorem, we obtain that 



-j(ip) dxdy = / / A-0 dxdy 



dtp 



dn dl + J ^^~ L ^y) d y 

-xp y (x,Y)dx- / xp x (0,y)dy, 

-L J-oc 

^-dl + cL, 
on 

where S~ is the top boundary of D~ , d/dn denotes the outward normal derivative 
at S~ , and dl means the line integration along S~ . The second equality uses the 
evenness and periodicity of ip in the re-variable and the last equality uses that 
ip y (x,y) — > — c as y — > —00. Since "f(ip) ^ we further obtain 

(5.23) - / ^-dl + cL^0. 



- dn 

On the other hand, since ip(x, y) — on S~ , it follows that 

g = (l + 7 ? 2 ( a ;))^ = -(l + 7 ? 2 ( a ;)) 1 /2|VV,| on S~ . 

Moreover, since (1 + ^(x)) 1 / 2 ^ r\ x (x) for — L < x ^ 0, the inequality (|5.23|) yields 
that 

rj x (x)\Vip(x, n(x))\dx < cL. 
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Finally, upon substituting of | Vt/>| by the Bernoulli's equation (|2.2b ) and integrating 
the above, we obtain 

{-2gn{-L)f/ 2 - (-2 OT (0)) 3 / 2 < cL. 
The assertion then follows by the use of the Bernoulli's equation again. □ 

Finally, the relative flow speeds at the crest and at the trough are bounded by 
the upstream speed of the underlying shear flow, and in turn, by the speed of wave 
propagation. 

Lemma 5.7 (Relative flow speed at the crest). Any nontrivial solution pair (A, w) € 
C's °f Ij2.13|) satisfies 

(\- 1/2 +w p (0,0))- 2 < X. 
Hence, the solution triple (c, r], tp) corresponding to(X, w) satisfies 

(5.24) ^(0,77(0)) < X = c 2 +2T 00 . 

Proof. Since w belongs to the nodal set M , it follows that w q = and w qq < on 
the half-line q = and — oo < p < 0. Recalling (|2. 12|) . subsequently, h q = and 
hqq < on q = and — oo < p < 0. Thus, (|2.9b ) along the half-line q = and 
— oo < p < reduces to the inequality h pp > — j(— p)h^. Further, 

~G?(b)) > - 2l[ ~ P) fo r-oo<p<0. 

Recalling that h p — > 1/c as p — > — oo, integration of the above inequality over 
(— oo, 0) then yields 

h~ 2 (0,0)-c 2 < 

The assertion then follows by the definition A = c 2 + 2roo. □ 

The same calculation carried out on the half-line q = ±L and p G (— oo, 0) leads 
to an analogous bound for the relative flow speed at the wave trough 

(5.25) V 2 (±i,77(±L)) > A = c 2 + 2r oo . 

5.2. Proof of the main results. Let C$ be the connected set in K + x C 3+a (R) x 
C 3+Q (ri^) of solution triples (c, rj, tp) of l]2.2[) . corresponding to the continuum C' s 
of solutions (A, w) of (|2.13|) via the transforms (|5.1|) and (|5.2p . 
By construction, C' s C C' s , if S > S'. Let 

<5>0 

and let C = sup 5>0 C$. 

For each S > 0, by virtue of Theorem 14.111 and Remark 14.131 at least one of the 
following holds: 

(1) there exists a sequence Wj)} C such that lim^oo Aj = oo; 

(2) there exists a sequence {(Xj,Wj)} C such that Xj are bounded but 
sup^- \uij I — > oo; 

(3) there exists a sequence {(Xj,Wj)} c such that Aj are bounded yet 
sup-^- \d p Wj | — > oo; 

(4) there exists a solution pair (A, w) £ C' s such that A + 2Fi n f = 5: 

(5) there exists a solution pair (A, w) £ C' s such that a _1 (A) + w p = 8 at some 
point in R; 
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(6) there exists a solution pair (A, w) G C' s such that w = 2X 4g 5 at some point 
on T. 

Proof of Theorem \2.1\ Due to (|5.ip . Alternative (1) implies that there exists a se- 
quence {{cj,r)j,ijjj)} C C$ for which lim^oo Cj = oo. This corresponds the first 
alternative in (ii) of Theorem 1 2. 11 

Our task is to give an interpretation of each alternative (2) through (6) in terms 
of the traveling speed c or the relative flow speed ip y to prove assertion (ii) of 
Theorem [2~T1 

Alternative (2). For each j and for each — oo < p ^ 0, by the nodal configuration 
of w it follows that d q Wj(q,p) > for — L < q < whereas by oddness of d q Wj it 
follows that w q (q,p) < for < q < L. Thus, w(q,p) attains in R its maximum 
somewhere along the line q — and its minimum somewhere along the line q = ±L. 
Therefore, this alternative implies that there exists a sequence {pj} in the interval 
(— oo, 0] such that 

either lim Wj(0,pj) = oo or lim Wj(±L,pj) = — oo. 

j— »oo j— >oo 

Suppose linij-Kx, Wj(0,pj) = oo. Since Wj(q,p) — > as p — > — oo for each j, 
furthermore, {pj} is bounded below. That is, there exists — oo < po < such that 

sup Wj(0,p) = sup Wj(0,p). 

On the other hand, for po ^ p ^ 0, 

supwj(0,p)= / d p Wj(0,p')dp' < |p |sup|<9pi<;j(0,j>)|. 

Thus, limj^oo sup p uij(0,p) = oo implies that linij^oo sup p d p Wj(0,p) = oo. That 
means, Alternative (3) must occur, as well. Below, we will show that in the case of 
Alternative (3), 

lim inf d y ij)j{x, y) = 

holds, where (cj,r]j,ipj) corresponds to (\j,Wj) via (|5.1[) and (|5.2p . The case 
limj^oo inf_ 00<p< o Wj(±L,p) = — oo is treated similarly. 

Alternative (3). Since Aj + 2r; n f > 5 and a _1 (Aj) + d p Wj > S for all j, it follows 
that a(p; Xj) > 5 1 / 2 for all j, and consequently, 

<5 < a~ l {\j) + dpWj < 8~ 1/2 + dpWj 

for all j, that is, 6 — S^ 1 ^ 2 < d p wj for all j. Therefore, limj_,oo sup^- d p Wj(q,p) = oo 
must hold. By (|5.2p and (|5.3p then it is readily seen that 

lim inf d y ipj(x, y) = 0, 
where (cj,r]j,ijjj) corresponds to (Xj,Wj) via (|5.ip and (|5.2p . 

Alternative (4). Let us choose a sequence {<5j} with £j — > 0+ as j — > oo and 
to each Sj choose (Xj,Wj) G C^. such that Aj + 2r; n f = Sj. We may assume that 
sup-^r dpWj < oo for all j ; otherwise, linij^oo inf n,, . d y ipj (x,y) =0 must hold by the 
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treatment for Alternative (3), where (c 3 , r\ 3 , ip 3 ) corresponds to (Xj,Wj) via (|5.ip 
and (J5T2J). 

Let us choose a sequence {pj}, — oo < p 3 ^ 0, such that 2T(pj) < 2r; n f + 5 3 for 
each j. It is straightforward that 

a(p j ;X j ) = (X j +2T(p j )) 1 / 2 <(2S j ) 1 / 2 , 

whence 

d p hj(q,pj) = o _1 (pj-; Aj-) + d p w 3 {q,p 3 ) ^ <5~ 1/2 + d p w 3 (q,p 3 ). 

Snce d p w 3 (q,p) is bounded, the right side increases unboundedly as j — > oo. Ac- 
cordingly, linij^oo supft d p hj(q,p) = oo, and 

lim inf d y ip 3 (x, y) = 0, 
where {cj,T]j,ipj) corresponds to (\ 3 ,w 3 ) via (|5.ip and (|5.2p . 

Alternative (5). Let us choose sequences {<5j} and {(Aj, tUj)}, {(9j,Pj)} such that 
Sj — » 0+ as j — > oo and (Aj, Wj) <G C' Sj , {q 3l p 3 ) € i? with 

dphj(qj,Pj) = a^ipf, \j) + dpWjiqj^j) = Sj. 

We may assume that p 3 is bounded below. Suppose on the contrary that pj — ► — oo 
as j — * oo. Since d p w 3 (q,p) — > as p — » — oo, it implies that 

and the right side tends to zero as j —>■ oo. It, in turn, implies \ 3 — > oo as j ' — ► oo. 
We may also assume that sup T \w 3 {q,p)\ is bounded. Otherwise, by the treatment 
of Alternative (2), lim J -_ >0O sup^ d y tpj{x, y) = holds. 

Let { (a,- , r)j , ;) } be the solution triples of (|2.2p corresponding to {(A^, w 3 )} via 
(|5.ip and (|5.2p . It is readily seen that d y ipj(xj,yj) — — oo as j — > oo, where !Ej = q 3 
and yj(x,pj) is the inverse of ij>j(x,yj). Under the assumption that po ^ Pj ^ 
for some — oo < po ^ 0, we claim that lim^oo sup po ^ p ^ d p h 3 {q,p) = oo, and 
correspondingly, 

lim sup dyipj (x, y) = 0. 

Suppose that d p h 3 (q,p) is bounded in po ^ P ^ for all j. Since 

/ d p h 3 (q,p')dp' = h 3 (q,0) - h 3 (q,p ) = iy(x) - y 3 {x,p ) < oo, 

and since rj 3 (x) is bounded for j, it follows that yj(x,po) is bounded, say ya < 
y 3 (x,po) for some yo. On the other hand, the pressure estimate (|5.14p yields that 

-\Vi(jj(x,y)\ 2 + gy -T(-ip j (x,y)) - -max(O,sup7(-0j))V , 3 < 

in . As j — > oo. However, inf^^^^) d y ipj(x,y) = -oo while all the other 
terms except for the first are bounded for j. A contradiction therefore proves the 
claim. 



Mi 
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Alternative (6). As is done for Alternative (5), let us choose sequences {Sj} and 
{Aj,u>} such that Sj — ► as j — ► oo and (Xj,Wj) G C' s , such that Xj — 2gvjj = 
somewhere on T. The nonlinear boundary condition i*2(A,u>) = then yields 

A + dpWjY = = — 

somewhere on T. As J — > oo, it implies that 

lim supd y ipj(x,y) = 0, 

where (cj,rjj,tpj) corresponds to (Xj,Wj) via (|5.1[) and (|5.2[) . 

In summary, there is a sequence of solution triples {(cj, rjj, ipj)} C C in the space 
R+ x C* 3+a (M) x C 3+Q (0,) such that 

either lim Cj = oo or lim sup cLV'j = 0. 

Vj 

This completes the proof. □ 

The remainder of this subsection is to refine the location of stagnation points in 
a "limiting" solution of (pETSJ) . 

Proof of Theorem \2.2\ If 7(7") ^ and "f'(r) for ^5 r < 00, then by virtue of 
Lemma [5.5[ the supremum of ipy m ^ is attained at the wave crest. Therefore, 
the only possible stagnation point is the wave crest. This completes the proof. □ 

In the finite-depth case [CS07, Var09], if the vorticity is non-positive (not nec- 
essarily monotone with depth), the only possible point of stagnation is shown to 
be the wave crest. In the infinite-depth case, the same result is expected. In other 
words, the monotonicity assumption in Theorem 12. 21 is expected to be removed. No 
proof is given presently, but some partial results are collected below. 

Since T(-ip) > for 7 (r) < 0, it follows by (f5TTT]) that 

|^(0,77(0))| < |VV(a?,tf)| for all (x,y) e fi,. 

Hence, if one can show that ip x = at the point of ip y = then one can obtain the 
desired result. Unfortunately, ip y or w p does not have a maximum principle, and 
hence it is not clear how to control the behavior of w p in terms of w q . 

Alternatively, if the vorticity is non-positive, Lemma l5. 31 states that ip y is mono- 
tone on the free surface. Thus, if one can show that the maximum of ip y in £l v 
occurs at the free surface then one can obtain the desired result. In the finite-depth 
case, it is shown [CS07,Var09] that if the vorticity is non-positive then ip y is mono- 
tone along the free surface, below the wave crest and the wave trough, and along 
the bottom. However, there is no sufficient control of the velocity at the infinite 
bottom. 

Proof of Theorem \2.3[ Assume the second alternative Iinx 7 _ ! . 00 sup^ d y tpj = oc- 

curs. We may choose a sequence {sj} and a sequence {(xj, yj)} such that Sj — ► 0— 
as j — > 00, (xj,yj) <E Sl Vj for each j and d y ipj(xj,yj) = Sj for each j. We may 
assume that {yj} is bounded from below; otherwise, the stagnation occurs at the 
infinite bottom, and Hindoo Cj = 0. 

Differentiating the Poisson equation Aipj = —j{ipj) in the y- variable leads that 

Adytpj + j'(ipj)dyxl)j =0 in to Vj 
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for each j. Let us introduce a sequence of functions 

(5.26) Wj(x, y) = d y ^(x, y) + s^"^ for (x, y) £ U, h , 

where (3 > is a constant such that (j 2 + j'(r) > for all ^ r < oo. It is 
straightforward that 

AWj + i{il)j)Wj = s.0 2 + i{^ ))e^ v - v ^ > in Q Vj 

and that Wj(xj,yj) = for each j. 

Since ^ 0, the weak maximum principle ensures that Wj in n Vj attains its 

maximum either on the surface or at the infinite bottom. On the other hand, 

Wj — > — Cj < as y — > — oo for each j. 

Therefore, Wj in Cl^. attains its maximum on the free surface y = f]j(x). Let 
—L ^ £j ^ be the maximum point of Wj in Cl n . . Since Wj(xj, yj) = 0, 
it follows that 

whence 

Since {y^} is bounded from below, by taking the limit as j — > oo we conclude that 

lim d y tpj(tj,Tij(Zj)) = 0. 

J—>OC 

That means ^ somewhere on the free surface becomes arbitrarily small. This 
proves (ii) of Theorem 12.31 

In case of a non-negative vorticity, ^ y (x,r](x)) is not necessarily monotone on 
the free surface, and one cannot expect that ip y = occurs at the wave crest. 
Nevertheless, ip y = cannot occur at the wave trough unless the free surface is flat. 
Indeed, if (x m , v( x m)) is the point of maximum horizontal velocity ip y on y = i](x) 
then by Bernoulli's equation, it follows that 

^(0,^(0)) + 2#7?(0) = ijjl(x m ,ri(x m )) +i/jy(x m ,ri(x m )) + 2gr)(x m )), 

whence 

tpl(x m , r](x m )) > 25(77(0) - rj(x m )). 

We assume, in addition, that i/j y (±L,T](±L)) is bounded along C. Then, by 
(|5.25p . the speed of wave propagation c is bounded along C, and the hrst alternative 
in (ii) does not occur. If, in addition, 7 is sufficiently small so that (|2.7[) holds, i.e., 

g + 7 (0)Mi> y (±L, V (±L))^0, 

then by Lemma 15.31 the relative flow speed ip y is monotone from crest to trough, 
and hence the second alternative in (ii) can be refined as 

lim d v ipj(a,ri(Q)) = 0. 

j— >oo 

This completes the proof. □ 



The smallness condition (|2.7p improves that (|5.20jl is non-positive [Hur06] , which 
involves solutions through the function n(x) in (|5.19|) . 
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6. Reformulation via a quasi-conformal mapping 

In the irrotational setting, Stokes [Sto80] proposed to use the (relative) veloc- 
ity potential (f>(x,y), defined in £l v as <p x — u — c, <p y = v, and the relative stream 
function tp(x, y) in order to study periodic traveling waves. By the conformafl hodo- 
graph transform (x, y) i— ► ((f), ip), the system (|2.2b )- (|2.2e ) is recast as Nekrasov's in- 
tegral equation [Nek51] (or Babenko's pseudo-differential equation [Bab87,BT03]), 
and an a priori bound for speed of wave propagation follows. 

With nontrivial vorticities, unfortunately, the velocity potential is not available. 
Nevertheless, under the no-stagnation assumption (|2.2a|l . the pseudo-velocity poten- 
tial and a quasi-conformal transform offer an alternative reformulation of (|2.2b V 
(|2.2b ). which share in common with the irrotational setting ( [LC25], for instance) 
some structural similarity. The development is adapted from [Zei73]. 

In preparation, let us rename the trivial solution (|2.10[) as 

c - utrty) = exp(r tr (V0) = (A + 2r(-^)) 1/2 , 
and let us make the ansatz 

u(x, y) — c = — exp(r tI . + r) cos 6, v(x, y) = exp(r tr + r) sin 9, 
or equivalently, 

2r _ ( U ~ C ) 2 + V 2 



e AT = —. — r^— , tanfl = . 

[Utr — cy U — c 

By construction, exp(2r tr + 2r) = (u — c) 2 + v 2 measures the kinetic energy density 
of the flow and 9 on the free surface y = r)(x) measures the angle that the wave 
profile makes with the positive horizontal direction. Note that exp(2rt r + 2r) > 
for a regular wave and a stagnation point corresponds to where exp(2rt r + It) = 0. 
In consideration of Stokes waves, r is required to be even and 2L-pcriodic in the 
cc-variable, and 9 is required to be odd and 2L-periodic in the x-variablc. 
Let us define the pseudo-velopcity potential 4>(x, y) in £1^ by 

(6.1) 4>x = W(x,y)^ y = W{u~c), ct> y = -W(x,y)ip x = Wv 

and </>(0, 77(0)) = for some function W(x,y). It is straightforward that (|2.2b ) 
dictates that the auxiliary function W satisfies 

(6.2) W x ^ x + Wyipy = Wity) in 0<y< V {x) 

It is reasonable to require that W(x, y) — ■> 1 as y — > —00. In addition, W is required 
to be positive, even and 2L-periodic in the x- variable. The complex potential (fi + iip 
is a p-analytic function. In the irrotational setting, W(x,y) = 1 everywhere in the 
fluid region and <j) + iip is & holomorphic function. Note that <f> is odd and ip is even 
in the a;- variable. Moreover, 

(f>{L, y) - (f>{-L, y) = / (u(x, y) - c)dx = -cL, 



independently of y. By the oddness of <f>, 

4>{±L, y) = TcL, 0(0, y) = for all < y < rj(x) 
and <j)(x + 2L, y) = 4>(x, y) — cL. 



1 Since —<j> is the harmonic conjugate of tp, the complext function (fi + itp is holomorphic in Q v . 
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Let us define the independent variables 
(6.3) q*=-<p(x,y) and p = -ip(x,y). 

They map the fluid region of one period {(x,y) G n n : — L < x < L} into the 
semi-infinite strip (—cL,cL) x (— oo,0) in the (q* , p) -plane and the free surface 
{(x,rj(x)) : —L < x < L} to the horizontal line segment (—cL,cL) x {0}. In what 
follows, let 

R* = {(q*,P) : -cL <q* <cL, -oo < p < 0} 

be the domain of one period in the transformed variables. 

The no-stagnation assumption, exp(2r tr + 2r) > throughout the fluid region 
guarantees the mapping (x, y) \— > (4>,ip) is quasi-conformal [Zci73]. In the irro- 
tational setting, the mapping is conformal. Furthermore, under this physically 
motivated stipulation, W is uniquely solvable in the C 1 class provided that ip is in 
the C 2 class. 

In view of r and 9 as functions of q* and p, under the change of variables 
(|6.3p . straightforward calculations yield that (|2.2b ) translates into the following 
inhomogeneous Cauchy-Riemann equations in the rectangle R as 

Wr q * + 6 P = 0, 

W9 q > +t p = + exp(-2(r tr + r)) 7 (-p). 

dp 

Furthermore, (|6.2p becomes 

e 2(Ttr+T) W 4 , = W-y(4>), W(q*,p) -> 1 as p -► -oo. 
Its unique solution is given explicitly in the (q* , p)-plane as 

(6.5) W(q*, P ) - cxp (^-y" P exp(-2(n r + r))7(-p)dp 

provided that exp(2r tr + 2r) > 0. The coefficient function W of (|6.4|) is a nonlocal 
operator involving the dependent variable r. In particular, (|6.4j) does not enjoy the 
maximum principle, and it may not be suitable for global existence theory. 
The change of variables in (|2.ip and (|6.ip is written in the concise form as 



W(u - c) Wv 



The back-transformation is given by 

x$ x^\ _ 1 (u — c —Wv 



,2/0 Vip) W{{u - c) 2 + v 2 ) \ v W(u-c) / 

A main advantage of this approach is that the nonlinear boundary condition on 
the free-surface takes a convenient form. Indeed, differentiation of the Bernoulli 
equation (|2.2d ) with respect to q*-variable yields the boundary condition 

r q . = . 9 A" 3/2 VF(0,r)" 1 e" 3r sin6l for p = 0. 

It uses that e Ttr ^ = A 1 / 2 . Since —Wr q * + 6 p = 0, furthermore, 

(6.6) P = g\- 3/2 e~ 3T sin6 for p = 0. 

The boundary condition (|6.6p is the same as that in the irrotational setting 
[LC25]. In the irrotational setting, 9 satisfies the Laplace equation A8 = in R* 
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and the boundary condition (|6.6p . By the sine series methods, then it leads to 
Nekrasov's integral equation 

sin |(s + 1) 



9( S ) 



1 

3^ 



log 



sin 2 (s — t) 



sm.e(t) 
v- 1 + f* sin 0(t')dt r 



-dt. 



where v = SgLc/nX 3 ^ 2 . In view of the positivity of the kernel, by multiplying the 
equation by sins and integrating over (0, tt) then yields that 

sin 9(t) sin t 



(s) sin sds = — 

6TT 



<- 



o v- 1 + J sin 9 (t')dt' 
9(t) sin tdt, 



-dt 



"3tt 

which implies that \i > 1/3tt. Therefore, an a priori bound for A (or v) is obtained 
( [BT03, Chapter 10], for instance). 
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